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ABSTRACT 

We develop a sophisticated model of FRB observations, accounting for the intrinsic cosmo- 
logical gas distribution and host galaxy contributions, and give the most detailed account 
yet of observational biases due to burst width, dispersion measure, and the exact telescope 
beamshape. Our results offer a significant increase in both accuracy and precision beyond 
those previously obtained. Using results from ASKAP and Parkes, we present our best-fit FRB 
population parameters in a companion paper. Here, we consider in detail the expected and 
fitted distributions in redshift, dispersion measure, and signal-to-noise. We estimate that the 
unlocalised ASKAP FRBs arise from z < 0.5, with between a third and a half within z < 0.1. 
Our predicted source-counts (“logN-logS”’) distribution confirms previous indications of a 
steepening index near the Parkes detection threshold of 1 Jy ms. We find no evidence for a 
minimum FRB energy, and rule out Emin > 10°85 erg at 90% C.L. Importantly, we find that 
above a certain DM, observational biases cause the Macquart (DM-z) relation to become 
inverted, implying that the highest-DM events detected in the unlocalised Parkes and ASKAP 
samples are unlikely to be the most distant. We do not expect our quantitative estimates in this 
region to be accurate until it is directly probed with localised FRBs. Since the cause of this 
effect is a well-understood observational bias however, it is guaranteed to be present to some 
degree. Works assuming a 1-1 DM-z relation may therefore derive erroneous results. 


Key words: radio continuum: transients — methods: statistical 


1 INTRODUCTION 


Fast radio bursts (FRBs) are radio transients of millisecond dura- 
tion and extragalactic origin (Lorimer et al. 2007; Thornton et al. 
2013). Their progenitors are unknown, with very many production 
mechanisms propoed (Platts et al. 2019). FRB surveys are providing 
increasingly large statistics with which to study the FRB population 
(Bhandari et al. 2018; Shannon et al. 2018; CHIME/FRB Collabora- 
tion et al. 2019a), including a handful of localised FRBs (Tendulkar 
et al. 2017; Bannister et al. 2019; Prochaska et al. 2019; Ravi et al. 
2019; Marcote et al. 2020). Furthermore, Macquart et al. (2020) 
have used localised FRBs as probes of the cosmological distribu- 
tion of ionised gas, illustrating their utility for cosmological studies 
(McQuinn 2014; Masui & Sigurdson 2015; Madhavacheril et al. 
2019; Caleb et al. 2019). Of the key questions surrounding FRBs, 
this work focuses on FRB population statistics. 

Studies of the FRB population are important both for under- 
standing the nature of FRBs themselves, and their use as cosmolog- 
ical probes. Typical fitted parameters include the FRB luminosity 
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function, e.g. minimum and maximum energies, and its shape; FRB 
spectral properties; and the source evolution. Studies including re- 
peating FRBs must also fit the distribution of repetition rates and 
allow for time-dependence between bursts from a single object. 


Linking FRB observations to the underlying true FRB popu- 
lation however is difficult. Connor (2019) review previous methods 
of studying the FRB population, and emphasise that accurate es- 
timates require accounting for the sensitivity effects of telescope 
beamshape, intrinsic burst width, and the dispersion measure dis- 
tribution p(DM|z) for a given redshift. In short, one must integrate 
over known or hypothesised intrinsic distributions of these variables, 
model observational biases, and then attempt to match observations. 
Doing so improperly will produce biased results. 


In what is usually seen as a different line of inquiry, cosmolog- 
ical studies using FRBs take advantage of their dispersion measure 
(DM), which integrates the column density of ionised gas along 
their line of sight. This encodes information on the diffuse gas in 
voids and galactic halos which is otherwise difficult to study. Mac- 
quart et al. (2020) have recently used the localised FRB population 
to constrain the total baryon density of the Universe and the degree 
of feedback. In that work, the authors analyse the probability dis- 
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tribution of observed dispersion measures, DM, given the redshift 
z of identified FRB host galaxies, p(z,DM). This controls some- 
what for the effects of the population of FRBs, which primarily 
affects the redshift distribution p(z). The authors do note however 
the potentially biasing effects of the FRB population as observed 
by FRB surveys, although a comprehensive treatement of such bi- 
asing effects is not performed due to the intrinsic error in using a 
small sample size (5-7). Many cosmological studies (such as helium 
reionisation) require many FRBs to both exist and be detectable at a 
redshift of z ~ 4 (Caleb et al. 2019), which is well beyond the most 
distant localized FRB to date at z ~ 0.6 (Law et al. 2020). 

Fundamentally, both lines of inquiry aim to study the intrinsic 
distribution of FRBs in z-DM space, p(DM, z). The only difference 
is the aspect of interest: population studies try to resolve a redshift 
distribution p(z) and treat the distribution p(DM|z) as a nuissance 
distribution, while cosmological studies aim to resolve p(DM|z) 
and attempt to remove the p(z) factor. They are thus fundamen- 
tally coupled problems. Unbiased estimates of the cosmological 
distribution of ionised gas require knowing the FRB population 
and the consequent biasing effects on measured dispersion mea- 
sures; and understanding the FRB population requires knowing the 
dispersion-measure distribution and its biasing effects on the mea- 
sured luminosity function. 

Caleb et al. (2016) provide the first comprehensive model of ob- 
servational biases on a simulated burst population, and FRBPOPPY 
(Gardenier et al. 2019) is being developed to include such effects. 
To date however, only Luo et al. (2020) have used this approach 
to fit population parameters. The authors study a sample of FRBs 
detected in the ~1 GHz band from Parkes, the Upgraded Molonglo 
Synthesis Telescope (UTMOST), the Australian Square Kilome- 
tre Array Pathfinder (ASKAP), Arecibo, and the Greenbank Tele- 
scope (GBT). The authors evaluate the validity of their model using 
MULTINEST (Feroz et al. 2009), which applies a Bayesian frame- 
work, and find a peak FRB luminosity L* of 29411210% ergs!, 


a differential power-law index of —1 T90 ae and a volumetric rate 


of 3791 104 Gpe~3 yr~! above 1042 ergs™!. 
In this work, we significantly improve upon FRB population 
models in the following manner: 


using an unbiased sample of FRBs from ASKAP and Parkes; 
using seven localised FRBs detected by ASKAP; 
correctly accounting for the full telescope beamshape; 

e using the measured signal-to-noise ratio in probability esti- 
mates; 

e including the intrinsic spread in the cosmological DM contri- 
bution due to large-scale structure and galaxy halos; 

e and allowing for redshift evolution of the FRB event rate per 
comoving volume. 


As with other population models, to make the problem tractable, 
we assume that the cosmological rate evolution, host galaxy DM 
contribution, burst width, and burst energy distributions are all in- 
dependent; and that FRB observations are random and uncorrelated, 
i.e. we do not model rapidly repeating FRBs. 

We begin by describing the ingredients to our model: a model 
of the DM distribution of FRBs as a function of redshift (Section 2), 
based on the model of Macquart et al. (2020); the intrinsic properties 
of FRBs, such as their burst width and the luminosity function, using 
a standard power-law description (Section 3); and the influence 
of observational effects such as beamshape and search sensitivity 
(Section 4). The method to combine these to calculate the expected 
“7-DM” distribution for observed FRB surveys is given in Section 5. 


In Section 6, we describe the data from ASKAP and Parkes to which 
we fit our model using maximum-likelihood methods. The best-fit 
FRB population parameters, and their uncertainties, are given in 
a companion paper (James et al. 2021). In this work, we present 
detailed comparisons to the observed DM, redshift, and signal-to- 
noise ratio distributions in Section 8, where we test for goodness- 
of-fit and search for deviations from expectations. Section 9 shows 
our estimates for the expected z-DM distribution of FRBs detected 
by ASKAP and Parkes. We summarize our results in Section 10. 
We attach in appendices a discussion of neglected effects in our 
modelling, and extra data for alternative source evolution scenarios. 


2 DISPERSION MEASURE DISTRIBUTION 


The distribution of dispersion measure, DM, of FRBs from a given 
redshift z, p(DM|z), is of both intrinsic interest, and is a nuissance 
factor in calculating the properties of the FRB population itself. 
Here, we use the method and parameters of Macquart et al. (2020). 
We model the DM of an FRB as 


DM = DMjsm + DMhalto + DMecosmic + DMhost: (1) 


with respective contributions from the Milky Way’s interstellar 
medium (ISM), it’s halo, the cosmological distribution of ionised 
gas, and the FRB host. In this work, we divide this into an ‘extra- 
galactic’ contribution, 


DMgG = DMeosmic + DMhost: (2) 
and a ‘local’ contribution, 
DMiocal = DMism + DMhato- (3) 


The ‘local’ contribution is subtracted from FRB observations, and 
thus all comparisons between expectations and measurements are 
made in terns of DMgg. This model slightly differs from that in 
Macquart et al. (2020), who model both DMhost and DMhalo using 
the same nuisance term, DMy. The distinction becomes important 
at large redshifts. 


2.1 DMism 


We use the NE2001 model (Cordes & Lazio 2002),! to estimate 
the Galactic contribution to dispersion measure. Since DM is an 
ingredient in the calculation of detection efficiency (see Section 4.2), 
the full integral for the FRB rate extends over the pointing direction 
as a function of Galactic coordinates, as discussed in Section A3 
and Eq. Al. Since most FRBs and FRB searches have been at 
high Galactic latitudes however, we use the mean value DM{sm to 
calculate the sensitivity for each survey, while using the individual 
values DMysm when calculating FRB likelihoods. 


2.2. DMhalo 


The exact contribution of the Milky Way halo to DM is uncertain, 
with estimates of order 10-80 pc cm~? (Prochaska & Zheng 2019; 
Keating & Pen 2020). FRBs have been observed down to a DM 
of little more than 110 pc cm~? (CHIME/FRB Collaboration et al. 
2019a) and 114 pe cm~? (Shannon et al. 2018), favouring the middle 
of this range and consistent with current estimations based on the 
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full set of observed DMs (Platts et al. 2020). We therefore use a 
value of DMhalo = 50 pecm™? in our default model. 

Deviations between our assumed values for DMysjyq and 
DMhalo Will be absorbed into our model for the host galaxy contri- 
bution. 


2.3 Cosmological DM 


We caution that symbols E, F, and a have different definitions in this 
section than in the remainder of this work. The notation regarding 
DMcosmic is derived from Macquart et al. (2020), and we preserve 
it for ease of reference to that work. 

The ‘cosmological’ contribution to DM, DMeogmic, can be 
understood as the DM incurred when an FRB is emitted at z at a 
random point in the Universe and propagates until the current epoch, 
z =0.A parameterisation based on detailed simulations (McQuinn 
2014) is given in Macquart et al. (2020), as a function of burst 
redshift z, as 


P(DMeosmiclZ) = DMcosmic) P(ADM)- (4) 


The expected value (DMeosmic) is calculated as per Ioka 
(2003); Inoue (2004): 


Ciie(z’)dz’ 
DMeosmic > 5 
$ ) o Ho(1+2')?E(z) m 


E(z) = YQm(1+2/)3 +Qq, (6) 


using the mean density of ions, ñe, and cosmological parameters 
relevant for the range 0 < z < 5: Ho = 67.4km s7! Mpc7!, and 
matter and dark energy densities Qm = 0.315 and Qa = 0.685 for a 
critical density Universe. See Macquart et al. (2020) and references 
contained therein for further details — cosmological parameters are 
taken from Planck Collaboration et al. (2018). 

The probability of deviations from the mean, p(Apm), is given 


by 
(Ape, - Co)* 


2:2 
2a ofM 


p(Apm) = AA exp (7) 


with a = 3, 8 = 3, and Cp being numerically tuned such that the 
expectation value of the distribution is unity. The degree of feedback 
F is reflected in the choice of opm = Fz~°->. In this work, we fix 
F = 0.32 based on results of (Macquart et al. 2020). 

The resulting distribution of DMgosmic,; P(DMeosmic!Z) is 
shown in Figure 1. 


2.4 DMhost 


The contribution of the FRB host galaxy (including the local en- 
virons of the FRB itself) to DM is highly uncertain. Some FRBs, 
most notably FRB 121102 and FRB 190608, show a large excess 
DM beyond what is expected from cosmological and MW contribu- 
tions, which cannot be explained by passage through an intervening 
galaxy along the line-of-sight (Spitler et al. 2014; Chatterjee et al. 
2017; Tendulkar et al. 2017; Hardy et al. 2017; Chittidi et al. 2020). 
Yet as noted in Section 2.2, many FRBs do not allow for a great 
deal of excess DM. Macquart et al. (2020) generically model this 
large spread using a log-normal distribution 


1 1 
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Figure 1. Distribution p(DM|z) of observation-independent DM param- 
eters, DMeosmic Only (top), and also including the best-fit distribution for 
DMhos derived in James et al. (2021) (bottom), as a function of redshift z, 
showing contours. 


We also correct the host contribution for redshift via 


DM; = DM‘ ost (9) 

ost 1 +z t 
In this work, we use Hhost ANd Chost as free parameters. Thus 
uncertainties in other DM contributions — including from our as- 
sumed value of feedback F — will be absorbed into these quantities. 


2.5 The intrinsic z-DM grid 


The probability distribution of observation-independent factors, 
DMeosmic + DMhost, is given in Figure 1. In this work, a linear 
grid in both z and DM space is used, with 1200 DM points spaced 
from 0-7000 in intervals of 5pc cm~’, and 500 in redshift from 
0.01-5. FRBs have their nominal local contributions, DMjoca] sub- 
tracted from their observed values of DM prior to evaluating their 
likelihood on this grid. In this model, only a small fraction of FRBs 
will have a DM very much larger than the mean. In particular, for 
z < 1.5, the majority of the spread in DM comes from the host 
galaxy, rather than the cosmological contribution. 
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3 FRB POPULATION 
3.1 Energetics 


Our model of the FRB population © is consistent with that used in 
the literature. We adopt a power-law distribution of burst energies 
E between Emin and Emax with integral slope y. We use ‘burst 
energy’ as the isotropic equivalent energy at 1.3 GHz, and assume 
an effective bandwidth of 1 GHz when converting between ‘per Hz’ 
and total quantities. 

In this model, the probability of a burst occurring above a 
threshold Ey, is a piecewise function, 


P(E>En) = 1 (E< Emin) 
p(E >En) = 0 (E> Emax) 


En Y _ Emax K 
Emin Emin 


p(E > En) = 


otherwise. (10) 


Observations show no evidence for a minimum FRB energy, and 
indeed the event rate is generally insensitive to Emin for y > —1.5 
(Macquart & Ekers 2018b). Thus we use a very low value of 
Emin = 1030 erg, which is several orders of magnitude below the 
minimum detected burst energy of all known FRBs — including 
SGR 1935+2154 at 1094-105 erg (The Chime/Frb Collaboration 
et al. 2020; Bochenek et al. 2020) — and treat this as a fixed pa- 
rameter. We re-examine this assumption in Section 8.3. 

Several other authors have used a Schechter function to model 
the FRB luminosity function (Lu & Piro 2019; Luo et al. 2020), 
which adds an exponential cut-off of the form exp(—E / Emax) to the 
luminosity function of Eq. 37. This is neither observationally nor 
theoretically motivated — the form of the Schechter function is used 
to model galaxy luminosities — and adds computational complex- 
ity, but avoids the unphysicality of a sharp cut-off. A preliminary 
investigation showed that using the Schechter function provided ef- 
fectively identical fits to the data. We therefore advocate for the 
simple power-law model out of simplicity. 

To calculate E, we convert between FRB energy E and ob- 
servable fluence F using 


47n D7 (z) 
E(F) = —— AF, (11) 
(1+ z)?2+e 
where a is the spectral index (F « v”), and Ay the bandwidth (here 
we use 1 GHz). Macquart et al. (2019) fit œ to 23 FRBs detected 
by ASKAP in Fly’s Eye mode, finding a = -1.5+0.2, Thus we use 
a default value of a = —1.5. We return to the interpretation of œ 
shortly. 


3.2 Population evolution 


The rate of FRBs per comoving volume will likely be a function 
redshift. While FRB host galaxies do not appear to be drawn from 
a population sampled proportionally to their star-forming activity 
(Safarzadeh et al. 2020), they certainly are not exclusively associated 
with very old galaxies in which star-forming activity has ceased 
(Bhandari et al. 2020a; Heintz et al. 2020). We therefore adopt 
the approach of Macquart & Ekers (2018b) and generically model 
the population evolution of FRBs as being to some power of the 
star-formation rate, i.e. 


®(z) 


®o (Fe (12) 


1 +z \SFR(0) 


5000 


Figure 2. Relative rate of FRB detections in z-DM space when ignoring 
beamshape and burst width effects for the best-fit FRB population parameters 
presented in James et al. (2021). All observational biases from Section 4 
have been ignored, and a constant detection threshold of 1 Jy ms is used. 
The Macquart relation, approximated by DM=(DM osmic )+€XP(Hhost)» 18 
also shown. The contours represent 50% (dotted), 90% (dash-dot), and 99% 
(dashed) of the FRB population. 


with Bg — and hence ®(z) — taking the units of bursts per proper 
time per comoving volume, i.e. bursts yr! Mpc”?. The factor of 
(1+z)7! converts between proper time in the emission and observer 
frames. We take SFR(z) from Madau & Dickinson (2014), 


1 2.7 
SFR(z) = ina e (13) 


1+ (+3) 
This model is useful in that n can be scaled as a smooth parameter. 
However it does not accurately model the source evolution should 
FRB progenitors originate from e.g. binary mergers with long delay 
times, as investigated by Cao et al. (2018). 
The total FRB rate in a given redshift interval dz and sky area 


dQ. will also be proportional to the total comoving volume dV, 


dv n (+z) D3) 
dad — “BE E@® 


(14) 


for angular diameter distance D 4, Hubble distance Dy = c/Hpo, 
and scale factor E (z) from Eq. 6. 

Applying this model with the best-fit FRB population param- 
eters derived in our accompanying work to the DM-distribution of 
Section 2 with a nominal threshold of 1 Jy ms produces the dis- 
tribution of FRBs shown in Figure 2. This ignores the important 
observational biases to be introduced in Section 4, and hence a 
quantitative analysis of the implications are left to Section 8. How- 
ever, it is clear that while at least 90% of FRBs will follow a 1—1 
DM-z relation (the Macquart relation Macquart et al. 2020), a sig- 
nificant minority will lie well above the DM-z curve. Indeed the 
highest DM events in a large sample are not likely to be the most 
distant. Consider DM> 3000 pc em™? in Figure 2. Such events can 
be produced by z ~3.2 lying on the Macquart relation. However, 
they must be the most intrinsically luminous FRBs to be detectable. 
At z ~1.6, observations probe a factor of ~ 4 further down the en- 
ergy distribution, allowing a greater number of events to be visible, 
and its high-DM tail may dominate the DM> 3000 pc cm? event 


MNRAS 000, 1-23 (2020) 


rate. This effect becomes more important for steeper luminosity 
distributions (large negative y) — this plot uses y = —1.2. 


3.3 Interpretation of a 


Many FRBs have a limited band occupancy (originally noted for 
FRB 121102 by Law et al. 2017), in which case the notion of 
a spectral index for an individual FRB has little meaning. In this 
case, the results of Macquart et al. (2019) can be interpreted as 
meaning either that there are more low-frequency FRBs, or that low- 
frequency FRBs are stronger. For an experiment with bandwidth 
similar to or less than that of the FRB bandwidth (which is the case 
with the data used here — see Section 6), the latter interpretation 
behaves identically to that of broadband bursts defined by a spectral 
index. However, the interpretation of an FRB population with a 
frequency-dependent rate does not. We denote this interpretation as 
the ‘rate interpretation’ of a, and that of Macquart et al. (2019) as 
the ‘spectral index’ interpretation. 

Under the spectral index interpretation of F(v) ~ v”, a neg- 
ative œ increases the detection threshold Fj, at high z due to the 
k-correction factor of (1 + z)~® through Eq. 11. This in turn de- 
creases the rate by a factor (1 + z)”® when Ey, << Emax through 
Eq. 10. Under the rate interpretation, the FRB population itself be- 
haves as ® ~ @(z,v) = ®(z)v®, and the k-correction therefore 
directly changes the rate, adding an additional factor of (1 + z)” 
to Eq. 12. Therefore, when y = —1, and Ey, << Emax, the two in- 
terpretations are identical. The situation becomes less simple near 
Emax, Which is frequency dependent under the spectral index inter- 
pretation, and constant under the rate interpretation — and the true 
behaviour may be more complicated than either result. 

Ultimately, we expect further observational data to be required 
to discriminate between the two scenarios, and consider both equally 
plausible for the time being. In this work, we present results using 
the spectral-index interpretation, but give additional data for the rate 
interpretation when constraining FRB population parameters. 


4 DETECTION THRESHOLD — OBSERVATIONAL 
BIASES 


FRB surveys usually calculate the fluence threshold above which 
FRBs would be detected using the radiometer equation, referenced 
to a 1 ms duration burst, using the sensitivity of the telescope at beam 
centre. This readily calculable value represents an unrealistic ideal. 
Bursts of longer duration will be harder to detect due to increased 
noise, while those viewed away from beam centre will be seen at 
less sensitivity. Furthermore, incoherent dedispersion searches will 
not perfectly match the shape of an FRB to the time-frequency 
resolution of the search, resulting in a lower detection efficiency. 

In this work, we model the effective fluence threshold Fj, as a 
function of nominal fluence threshold at 1 ms F4, beam sensitivity 
B (normalized to a maximum of 1), and an efficiency factor due to 
burst duration, 77, as 

Fh = (15) 

This results in a theoretical distribution of bursts in z-DM space, 
p(z, DM), as 


PDM) = f dB f anpi DMIF tn B08) pC. (16) 


where p(z, DM|F;,) is the distribution at a fixed threshold (Fig- 
ure 2), Q(B) is the region of sky at which the beam sensitivity is B, 
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and p(7) is the probability that burst properties lead to a total de- 
tection efficiency 7. The effects of these two factors are investigated 
in the following sections. 


4.1 Beamshape 


A telescope’s beamshape is usually represented as the relative sen- 
sitivity B as a function of the sky position Q relative to boresight, 
B(Q), such that B(O) = 1. The beamshape is often approximated 
as a Gaussian or Airy function, although precise measurements of 
B can become important when attempting to localise FRBs de- 
tected in multiple beams, or estimating the relative rate of single- 
vs multiple-beam detections (Vedantham et al. 2016; Macquart & 
Ekers 2018b). 

For the purpose of estimating the number of FRBs detected 
however, the ‘inverse beamshape’, Q(B), which describes the 
amount of sky Q viewed at any given sensitivity B, becomes more 
relevant (James et al. 2019a). Most calculations of FRB rates have 
characterised a telescope beam as viewing out to the FWHM at 
full sensitivity, i.e. Q(B) = Qpwym6(B - 1) (e.g. Thornton et al. 
2013; Bhandari et al. 2018). Others have used a Gaussian approx- 
imation for the beamshape (e.g. Lawrence et al. 2017), which is 
equivalent to Q(B)dlog B = const. We here analyse the suffi- 
ciency of these approximations, using for the Gaussian approxi- 
mation o = (FWHM/2)(2 log 2)9-5, where the full width at half 
maximum (FWHM) assumes an Airy disk, i.e. HPBW=1.22A/D 
for wavelength at central frequency A and dish diameter D. 

ASKAP FRB observations have varied the observation fre- 
quency and configuration of beams formed from ASKAP’s phased 
array feeds (PAFs). However, the majority of both fly’s eye and 
incoherent sum observations have used the ‘closepack36’ config- 
uration at a central frequency of 1.296 GHz. We therefore use the 
beamshape derived in James et al. (2019a). In the case of the Parkes 
multibeam, we use a central frequency of 1.382 GHz, and the sim- 
ulations of K. Bannister (published as Vedantham et al. 2016) and 
A. Dunning (referenced as ‘private communication’ by Macquart 
& Ekers 2018a)), which produce equivalent results for Q(B). This 
also allows us to conclude that while the ability to localise FRBs 
detected in multiple beams may be limited by systematic uncer- 
tainties in the beamshape (Macquart & Ekers 2018a), the inverse 
beamshape B(Q) is robust against such certainties, since it does not 
care about where on the sky any given patch of sensitivity is located. 

Figure 3 shows the resulting ‘inverse beamshape’ function 
Q(B). This is compared to the equivalent Q(B) when using the 
Gaussian and Qfwym approximations. Since implementing the full 
function Q(B) in the calculation of the z-DM distribution is nu- 
merically expensive, we investigate the accuracy of reducing Q(B) 
to a small number of values. A set of such values is also shown 
in Figure 3. The accuracy of all approximations is assessed against 
the full beamshape, by comparing the total predicted number of 
events and the mean value of DM to those calculated for the full 
beamshape, and also assessing the maximum difference between the 
p(DM) curves, as shown in Figure 4. Table 1 lists the resulting er- 
rors. This is evaluated for the best-fit set of parameters found in our 
accompanying paper — however, a brief investigation has shown 
that results are not sensitive to the assumed parameters within a 
reasonable range. 

We find that using five values of B for ASKAP, and ten for 
Parkes, achieves an approximate DM distribution with an error in 
p(DM) of less than 1%, and an error in the mean value DM of 
only 0.2%. In contrast, using the FWHM approximation, which is 
the standard in the current literature, results in O ~10% deviations 
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Figure 3. Inverse beamshape Q(B) as a function of beam sensitivity B for 
ASKAP (closepack36 configuration at 1.296 GHz) and the Parkes Multi- 
beam (1.382 GHz). Shown are the best measurements (lines), FWHM ap- 
proximation (crosses), the numerical approximations used in this work (cir- 
cles), and the Gaussian beamshape (flat horizontal lines). The points for 
the FWHM and numerical approximations represent 6 functions — in the 
FWHM case, the true amplitude is shown, while in the numerical approx- 
imations, the amplitudes are renormalised by their spacing in log B for 
comparability with the best measurements. 
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Figure 4. Top: dependence of expected FRB DM distributions (pc cm~?) on 
beamshape models for three different FRB surveys. The beamshape models 
considered are shown in Figure 3: the‘ FWHM’ approximation, the ‘Full’ 
beamshape Q(B), a numerical approximation used in ‘this work’, and a 
Gaussian beamshape. All curves for each telescope are normalised such 
that the distribution for the full beamshapes peaks at unity. Bottom: the 
difference between the ‘Full’ beamshape and that found when using the 
‘FWHM’, ‘Approx.’, and Gaussian approximations. 


in the DM distribution and its mean value, and pushes the mean 
towards higher values. The total expected detection rate found when 
using the FWHM approximation is almost double that when found 
when using the full ASKAP beam, but there is only a 20% excess 
for Parkes. When assuming Gaussian beams, a huge excess in the 
total rate of ASKAP bursts is predicted, since this does not account 
for the closely packed, and thus overlapping, beams. We note that 


Table 1. Percentage errors in the total FRB rate and mean DM value, DM, 
and maximum deviation |6DM|max, when using different beam approxima- 
tions when compared to that found for the full beamshape function Q(B), 
for each of three FRB surveys considered. 


Survey Approximation Rate DM = |6DM|max 
FWHM +76 +6 0.1 
ASKAP/FE This work +17 +0.2 0.009 
Gauss +513 -3 0.07 
FWHM +73 +8 0.11 
ASKAP/ICS This work +5.8 +0.2 0.005 
Gauss +440 -4 0.08 
FWHM +20 +14 0.15 
Parkes/Mb This work +2.8 +0.2 0.005 
Gauss +160 +1 0.011 


uncertainties in the true ASKAP beamshape due to the calibration 
procedure (see James et al. 2019a) are less than the errors introduced 
by our numerical approximation. 

In the case of Parkes/Mb, the excess rate when using a Gaus- 
sian beam is due to outer beams being less sensitive than the central 
beam at which the sensitivity is usually calculated. However, the 
Gaussian beam approximation accurately estimates DM and the 
shape of p(DM) distribution. This suggests that even very complex 
beamshapes, such as that of the Canadian Hydrogen Intensity Map- 
ping Experiment (CHIME), could be included in our model in a 
relatively simplified manner. 


4.2 Detection efficiency 


We model the threshold at which an FRB of fixed fluence F can be 
detected as scaling with the square root of its effective width, weff, 
relative to the nominal width of 1 ms, using an efficiency factor 7: 


n = = (17) 
Fin(Werr) 
_ Werf \~9-5 
s(a (18) 


The effective width is modelled as per Cordes & McLaughlin 
(2003), being a function of its intrinsic duration wint, scattered 
width Wscat, DM smearing within each frequency channel wpm, 
and the time-resolution of the search Wsamp: 


= 2 2 2 2 
Wef = ATA + Weoat + WDM + Wsamp- (19) 


Often, the scattered width and intrinsic width are indistinguishable, 
and their separation only becomes important for telescopes observ- 
ing at different frequencies. We therefore define the ‘incident’ width, 
Wines aS 


Wine = Wint + Wscar (20) 
An alternative model is presented by Arcus et al. (2021), which is 
based on fits to simulated ASKAP and Parkes FRBs. Since it is 
not clear how the fit parameters translate to the general case, and 
because we wish to present a broadly scalable model, we do not use 
their formulation. We remark rather that the widely used model of 
Eq. 19 should be investigated further. 

In order to model the distribution of wince, P(Winc), We use a 
log-normal distribution, 


1 (= Winc—l08 Hw \ 
7 M 


Wine log ow 


— e 21 
log cw (27)9-5 on 


P(Winc)4Wine = 
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Figure 5. The effect of intrinsic burst width on survey sensitivity. Blue: in- 
trinsic distribution of burst widths wint (peak set to unity). Orange/green/red: 
the simulated width distributions expected from the three FRB surveys after 
accounting for observational bias. Two normalisations are shown: the upper 
curves have been normalised to a peak probability of unity, for comparison 
with the fit to observed values by Arcus et al. (2021) (shown in black); 
while the lower curves are normalised relative to the rate for w = 0, and 
thus illustrate the relative reduction in event rate as a function of intrinsic 
width. The upper curves for ASKAP/FE and ASKAP/ICS coincide almost 
identically with the results of Arcus et al. (2021). 


Table 2. The effects of different assumptions on the intrinsic burst width 
distribution at 1.3 GHz for three FRB surveys. Given are the total detection 
rate (normalised to a rate of unity at w = 0), and mean DM DM. The width 
distributions are parameterised via Eq. 21, with given values of uy, and 
Ow corresponding to no intrinsic width, the observed width distribution 
of Arcus et al. (2021), when accounting for observational bias to find the 
intrinsic distribution, and when numerically approximating (~) the intrinsic 
distribution with a few characteristic values. 


Parameter Hw Ow ASKAP/FE ASĶAP/ICS Parkes 
0 0 1 1 1 

Rate 207 207 0.46 0.51 0.20 

ae 549 246 0.27 0.30 0.11 
~5.49  ~2.46 0.27 0.30 0.11 

0 0 263 371 488 

pm 207 207 286 397 724 
5.49 2.46 293 401 726 

~5.49  ~2.46 292 400 724 


We do not include any DM or z dependence in the width distribution 
— see Appendix A4 for further discussion on this topic. 

Unlike Luo et al. (2020), we do not include fits of the model 
parameters Wy and oy as part of our general fitting process. Rather, 
we use the low correlation between uw, oy and other parameters to 
first fit for these values using a preliminary parameter set, and then 
check that the fit is still valid for the final parameter set presented 
in our companion paper (James et al. 2021). 

Arcus et al. (2021) use Eq. 21 to model the observed width 
distribution p(winc) of ASKAP and Parkes FRBs, finding uw = 
2.67 ms and oy = 2.07. We instead use Eq. 21 to model the intrin- 
sic width distribution, and vary fy, oy until the simulated width 
distribution of the ASKAP/FE survey matches the parameterisation 
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Figure 6. Effect of using different distributions of the intrinsic width wint on 
the expected DM distribution of FRBs, p(DM), from the FRB surveys con- 
sidered here: ASKAP/FE (left-most), ASKAP/ICS (centre), and Parkes/Mb 
(right-most). Assuming no intrinsic width (blue solid), using a log-normal 
distribution of widths with parameters from Arcus et al. (2021) (green dot- 
ted), using a log-normal distribution with parameters derived in this work 
(red dashed), and numerically approximating the latter (black dot-dashed). 
The distributions are normalised to a peak of unity for illustrative purposes. 


of observed widths by Arcus et al. (2021). We obtain uw = 5.49 
and oy = 2.46, and then proceed to use these values in further 
calculations to optimise FRB population parameters. This is shown 
in Figure 5, with total rates, and mean estimated DM values, given 
in Table 2. Finally, we re-evaluate the fits using the final optimal 
set of parameters we present in our companion paper (James et al. 
2021). 

Comparing the intrinsic (blue) and observed (black, coloured) 
distributions in Figure 5, modelling the intrinsic FRB rate, and 
accounting for observational bias, correctly reproduces the observed 
FRB width distribution as estimated by Arcus et al. (2021). The 
effect of observational bias is clearly reflected in the total expected 
FRB rate, with high-width bursts much less likely to be detected. 
The magnitude of this effect, shown in Table 2, depends on the 
significance of the other terms in Eq. 19. When the time resolution 
is poor (Wsamp large), the effect of a large intrinsic width is less — 
thus the reduction in rate for ASKAP is less than that for Parkes 
surveys. A second-order effect is that more-sensitive surveys, which 
probe further into the Universe and see FRBs with on-average higher 
DMs, are also less sensitive to wint, since Wpm is larger. Hence 
the reduction in rate for ASKAP/ICS is slightly lower than for 
ASKAP/FE, while the greatest effect is for the Parkes/Mb survey, 
where the intrinsic FRB width reduces the number of detected FRBs 
by a factor of 10. 

A consequence of this is that the true number of high-width 
FRBs will be very difficult to estimate, so that our log-normal 
model is effectively untested beyond w = 10ms. Thus while we 
estimate that ASKAP/FE and ASKAP/ICS miss ~ 70% of FRBs 
due to their intrinsic width, and Parkes/Mb 90%, we do not consider 
this quantitatively reliable — there may be virtually any number of 
high-width events remaining to be detected. 

Whatever the lost rate, losses will preferentially arise from the 
nearby Universe where wpy is low. Including the width distribution 
therefore increases the mean expected DM, DM. This effect is small 
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(~ 10%) for ASKAP observations, but more significant for Parkes 
(~ 30%). 

Finally, we note that while including the width distribution is 
clearly very important, the details matter less. Using the parameters 
of Arcus et al. (2021), or approximating the true distribution with 
a small number of points for computational efficiency, produces 
an almost identical value of DM. This also means that the loss of 
efficiency to high-width bursts for the HEIMDALL? search pipeline 
found by Gupta et al. (2021) — which is commonly used in Parkes 
FRB searches — is insignificant to our modelling. 


4.3 Numerical implementation 


The integrals over B and 7 in Eq. 16 are numerically expensive. Fur- 
thermore, we have shown above that approximating the beamshape 
with 5—10 values, and the width distribution with five, allows for a 
very good approximation to the continuous distributions. Therefore, 
we approximate these continuous distributions with these discrete 
distributions, i.e. 


NB 
Q(B) ~ X 916(B - Bi) (22) 


Nw 
p(n(DM,w)) ~ >) pid(n(DM, ww) — (DM, wi)),23) 


recalling that 7 is a function of both DM and w through Eq. 18 and 
19. Therefore, the z-DM distribution of Eq. 16 becomes a weighted 
sum over individual distributions at fixed observational thresholds, 


Ng Nw 


p(z,DM) =" >) p (& DM, Fin(Bi.7(DM, w;))) Qipj- (24) 
i=l j=l 


5 METHODOLOGY 


The ingredients described above are implemented in PyTHon. Here, 
we describe its effective implementation, which is identical in ap- 
proach to the method proposed by Connor (2019) and implemented 
by Luo et al. (2020), even if the exact details differ. 

We consider FRB data from some number of independent 
FRB surveys. The total likelihood £ of the outcome of multiple 
independent FRB surveys is simply a product of their individual 
likelihoods, £j, 


Nourveys 
£= || £ (25) 
i 
for surveys i = 1 . . . Nsurveys. We describe L; as the product of three 
independent terms: 


Ni 
Li = pn (Ni) | | Pamz(DMj.z;)ps(sjIDMj.z;). 20 
j=l 


Here, pn is the probability of detecting N; FRBs in survey i, while 
Pamz İs the likelihood of the j th FRB from survey i being detected 
with dispersion measure DM; and, when applicable, at redshift z j. 
We also include the probability ps of an FRB being detected with 
fluence F a factor s above the fluence threshold Fip given it was 


2 https://sourceforge.net/p/heimdall-astro/wiki/Home/ 


observed with properties DM ;, zj. Each of these terms is described 
independently below. 

It is also possible to separate further terms in Eq. 26. As noted 
by Vedantham et al. (2016), for a multibeam instrument, the relative 
likelihood of a single- vs multi-beam detection, and the relative 
likelihood of detection in different beams of varying sensitivity, are 
functions of the FRB fluence distribution. Such measures are only 
relevant when the true FRB fluence F cannot be resolved, as is 
common with FRB detections by the Parkes multibeam. When F 
can be reconstructed, as is the case for CRAFT FRB detections with 
ASKAP (Shannon et al. 2018), then the survey acceptance to that 
particular FRB can be calculated exactly, and p(s) in Eq. 26 can be 
written in terms of p( F). Similarly, we will not include the observed 
value of an FRB’s width when evaluating pgm, and ps, with only 
the overall width distribution being accounted for. We consider that 
adding such terms will yield only a small increase in analytic power 
for resolving the FRB population, at the cost of a large increase in 
complexity. Thus they are ignored — however we do acknowledge 
that we are discarding a small amount of information by doing so. 


5.1 Probability of N detections, p,(N) 


Ignoring the correlations caused by repeating FRBs (see the discus- 
sion in Section 5.4), the total number of observed FRBs in survey 
i, N; comes from a Poisson distribution, 


(Ni)? exp(— (Ni)) 

Ni! i 
where (N;) is the expectation value of N;. The calculation of (N;) 
is the heart of the problem that we address in this work, since it 
must necessarily incorporate all relevant properties that affect the 
detection rate. 


Combining the dependencies in Sections 2-4, (N;) is calcu- 
lated as 


(Ni) = TiRi, (28) 


dV 
J dz D(z) T 


Pn(Ni) (27) 


Ri 


dDMp(DM|z) (29) 


Il 


J aoB) f dwp) pE > E(B, w2 DM). 


Here, T; is the survey duration, which multiplies a rate R; to pro- 
duce the total expected number of bursts. ®(z) is the FRB source 
evolution function (Eq. 12), dV(z)/dQ/dz is the comoving volume 
per steradian per redshift interval from Eq. 14, p(DM|z) is the ex- 
tragalactic DM distribution found by convolving Eq. 4 and 8 (shown 
in Figure 1); Q(B) is the beamshape discussed in Section 4.1 and 
approximated as per Eq. 22; p(w) is the width distribution of Eq. 21, 
discretised as per Eq. 23; and p(E > E) is the cumulative energy 
function of Eq. 10. The dependency of this threshold E on the 
parameters (B,w,z,DM) is encapsulated in Eq. 11, Eq. 15, and 
Eq. 18. 

For some surveys, no controlled survey time 7; is available, 
and this term is simply set to unity in Eq. 26. However, R; can 
be calculated regardless of knowledge of 7;. For those surveys with 
known 7;, the most likely value of the lead constant in the population 
function of Eq. 13, ®p, can be estimated without recalculating the 
integral of Eq. 29. 


5.2 Calculating pamz 


The probability of an FRB being observed with a given dispersion 
measure DM and redshift z is given by the appropriate integrand of 
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Eq. 29, 


dV(z) 
dQdz 


pane = Rao pom) dBQ(B) (30) 


f wooo > Em(B, w, z,DM)). (31) 


For FRBs with no measured host redshift, the relevant quantity is 


Pdm = f Pamz@z, (32) 


and pgm replaces pqmz in Eq. 26. The rate R; is used as a normalising 
factor in Eq. 31, so that 


J J PamzdDMdz = J PamdDM =1. (83) 


The shape of pgmz in z-DM space is a primary quantity of interest 
in this work. 


5.3 Calculating p,(s) 


The measured fluence F of an FRB also holds information on the 
FRB population. However, in many telescope systems — and no- 
tably for Parkes (Macquart & Ekers 2018a) — F is not directly 
measured, since the location of the FRB in the beam is not known. 
Furthermore, we are interested in p(F'| Fy), i.e. the probability of 
measuring F given an FRB has been detected at threshold Fip, 
which itself has complex dependency through Eq. 15. 

This difficulty can be overcome by noting that the signal-to- 
noise ratio, SNR, is a readily observable parameter for an FRB, 
and most FRB-hunting systems use a well-defined threshold SNR, 
SNR, to distinguish FRBs from noise. As per James et al. (2019b), 
who base their work on Crawford et al. (1970), we define 


SNR 
= 34 
SNRin Bo 
F 
= = (35) 
Fih 


where Fin is the fluence threshold to that FRB. As detailed in Sec- 
tion 4, F is a function of the burst DM, width, and the location in 
which it is observed by the telescope’s beam, so that neither term in 
Eq. 35 is known. However, the ratio is preserved in the measurable 
quantities of Eq. 34, making s a very useful observable. 

The probability p(s) of observing s in the range s to s + ds 
given that an FRB has already been observed is 


2 a 
Ps(s|DM,z) = — | waw) f awp)" j 
= 8f dV (z) 
Ps(s|DM) = — feo Tad; p(DM|z)dz 


J eae f awpow Je aper in) (36) 


for localised and unlocalised FRBs respectively. In the integrands, 
p(sFp) is the probability of detecting a fluence F = sFy, in an 
interval between s and s + ds, where Fy, depends on B, w, and DM 
as per Eq. 15. 

The probability p(sFy,) can be found from Eq. 10. Given that 
such an FRB has been observed at all, the integral distribution of 
E given that an FRB has been detected can be found by replacing 
Ep with E and Emin with E. Differentiating by E produces the 
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probability amplitude 


g \r 
dp (EobslEth) 7 Y (£) 
dE Eh 1- (3) 
En 
y oe 


= a = l Epa) (37) 


The value of Ep can be found as a function of z by inserting Fih 
into Eq. 11; thus p(sFy) is equivalent to p(sEm(z, Fh)). Relating 
dp/dE from Eq. 37 to the required dp/ds of Eq. 36 produces 


dp (Eopbs|Etn) dp(Eobs|Etn) dE dF 


E 38 
ds dE dF ds (38) 
dp(Eops|Etn) E 
p( obs th) En Fa (39) 
th 
sy] 
at es 4 
J Ea” (40) 
En 


Inserting Eq. 40 into Eq. 36, and integrating over the appro- 
priate dimensions, produces the desired ps. 


5.4 What about repeating FRBs? 


By writing the individual burst probabilities as being independent 
in Eq. 26, and assuming that the number of detected bursts follows 
a Poisson distribution in Eq. 27, we ignore the potential of FRBs to 
repeat. While the fraction of the FRB population which is observed 
to repeat is a current topic of debate, it is certain that many do. 
Formally, the FRB population described in Section 3 represents all 
bursts, rather than all FRB emitters, and the summation of Eq. 26 
runs over all detected bursts. The distinction becomes irrelevant 
for distant, rarely repeating sources for which only ever zero or one 
bursts will be observed. For bright, nearby repeaters, the probability 
of having such an object in a survey’s field of view is rare, especially 
when burstiness and/or periodicity is included (Oppermann et al. 
2018; Chime/Frb Collaboration et al. 2020; Rajwade et al. 2020; 
Cruces et al. 2020), and an observation of zero bursts will be more 
likely than that estimated by Eq. 27. Conversely, the probability of 
observing many bursts will also be high, with observations of single 
bursts being much rarer than otherwise expected. 

We note that the only FRB survey to ever observe an FRB repeat 
in an unbiased way are the observations by CHIME (CHIME/FRB 
Collaboration et al. 2019b; Fonseca et al. 2020a,b) — all other 
repeating FRBs have been discovered in targeted follow-up obser- 
vations. This suggests that the majority of FRB observations can 
safely be classified as being in a ‘one burst per progenitor’ regime, 
regardless of the true fraction which are actually repeating objects. 
For these, our approach should be valid. We revisit this assumption 
in Section 8.4. 


6 SURVEYS 


Estimates of the FRB population have been made using data from 
many telescopes, which are often drawn from FRBCAT (Petroff 
et al. 2016). Due to the large number of FRBs they have detected 
and published, results from Parkes and ASKAP remain the most im- 
portant, and we focus on these instruments here. Other important in- 
struments we wish to examine in future works include the Upgraded 
Molonglo Observatory Synthesis Telescope (UTMOST), and the 
Canadian Hydrogen Intensity Mapping Experiment (CHIME). 
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The sensitivity of an FRB survey — and hence the functions 
PN, Pzam, and ps from Section 5 — depends on the local contri- 
bution to DM, and hence varies with the value of DMysq. Since 
this fluctuates pointing-by-pointing, in theory these functions must 
be recalculated for every single pointing direction, which becomes 
computationally prohibitive. Evaluating pzqm and ps however for 
the measured DMgg and s of an FRB is much quicker. This moti- 
vates grouping FRB observations not just by telescope, but also by 
other observational properties, such as Galactic latitude. Here, we 
use five groups, as described below. 


6.1 Parkes 


All Parkes FRBs published so-far have used the multibeam (‘Mb’) 
receiver (Staveley-Smith et al. 1996). However, a new ultra- 
wideband receiver is now in place (Hobbs et al. 2020), which is 
being used for FRB searches and follow-up observations. We there- 
fore refer to results from Parkes as “Parkes/Mb” to distinguish this 
from future works. 

Of the many Parkes FRB discoveries, we consider only those 
by the High Time Resolution Universe (HTRU; Keith et al. 2010; 
Thornton et al. 2013; Champion et al. 2016; Petroff et al. 2014) 
and Survey for Pulsars and Extragalactic Radio Bursts (SUPERB; 
Keane et al. 2017; Bhandari et al. 2018) collaborations to have an 
unbiased estimate of observation time, T;. This is because their 
observation time and pointing directions were pre-determined, and 
the results published regardless of outcome. Other results suffer 
from publication bias whereby non-detections are less likely to be 
published. Thus, while their discovery can contribute to individual 
FRB likelihoods via pqmz and ps, no well-defined observation time 
T; exists for use in Eq. 29, and the term p;(N;) must be neglected 
in Eq. 26. Thus they cannot contribute to estimates of the total FRB 
rate. 

Since the distribution of DMysjy affects telescope sensitiv- 
ity, surveys at low Galactic latitudes have significantly reduced 
sensitivity compared to those at high latitudes, and the full dis- 
tribution of time spent at each DMysyjy must be accounted for. In 
particular, DMysyy will vary significantly for each pointing at low 
latitudes, making estimates numerically taxing. We therefore in- 
clude only FRBs detected at mid (19.5° < |b| < 42°) and high 
(42°le|b| < 90°) Galactic latitudes. This criteria leaves 12 FRBs 
detected in a total of 164.4 days by HTRU and SUPERB (Bhandari 
et al. 2018), and another 8 FRBs by other groups with no reliable 
observation time. A full list is given in Table 4. 

Early searches for FRBs with Parkes used a sparse grid of DMs 
and arrival times, resulting in sensitivities that would fluctuate by 
+15% (Keane & Petroff 2015b). This was corrected with the use 
of HEIMDALL’, which has been used to (re)process the data from 
HTRU and SUPERB. While early HTRU searches extended only to 
DM=2000 pc cm~? (Thornton et al. 2013), latter searches extended 
this to 5000 pc cm3; and while the SUPERB ‘F’ pipeline looks for 
FRBs with DM< 2000 pc cm, the SUPERB ‘T’ pipeline extends 
the search to 10,000 pc cm~’. Thus we treat all Parkes FRB searches 
as fully covering DM space. 

For the Parkes multibeam, nominal sensitivity to a burst at 
beam centre is F, = 0.5 Jy ms to a 1 ms duration burst (Keane et al. 
2017) — this is an approximation, since different FRB searches 
used slightly different values of SNR. We neglect the effect of 
1-bit sampling with early searches for FRBs with Parkes, which 
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would have slightly degraded the sensitivity of these observations. 
These and other properties are summarised in Table 4. 


6.2 ASKAP 


The Commensal Real-time ASKAP Fast Transients (CRAFT) group 
have performed several FRB surveys with ASKAP. The majority 
of ASKAP FRBs have been observed in single-antenna (“Flye’s 
Eye”, or “FE”) mode during the ‘lat50’ survey, i.e. while observing 
Galactic latitudes of |b| = 50° + 5° (Bannister et al. 2017; Shannon 
et al. 2018). Twenty FRBs have been initially reported (Macquart 
et al. 2019; Bhandari et al. 2019), with a total recorded data time 
of 1274.6 antenna-days duration (James et al. 2019a). A further six 
FRBs have been detected in a variety of surveys (Macquart et al. 
2019; Bhandari et al. 2019; Qiu et al. 2019), of which four satisfy 
the Galactic latitude requirement. These are listed in Table 5. We err 
on the side of caution and do not assume a known observation time 
for this last category, since several other FRB searches outside the 
lat50 survey have been performed and were not reported. As with 
Parkes, p;(N;) in Eq. 26 is only evaluated for the former category. 

All ASKAP/FE searches have used the same frequency range 
and time/spectral resolutions, as given in Table 3. The beamshape 
and threshold for this survey are given in James et al. (2019a). 

ASKAP has recently been observing in incoherent sum mode 
(ICS), with voltage buffers used in offline analysis to localise FRBs 
to sub-arcsecond precision (Bannister et al. 2019). Follow-up ob- 
servations with radio and optical instruments have determined the 
redshifts of the host galaxies of each FRB, allowing the DM-z grid 
to be directly probed for the first time. This mode has undergone an 
extended period of commissioning, with the number of telescopes, 
observation frequency, and time resolution of the search all varying. 
The total observation time is difficult to estimate, again precluding 
the use of observation time in this survey’s likelihood calculation. A 
comprehensive analysis would involve recalculating the z-DM grid 
for each and every observed burst. For reasons of computational 
efficiency, we instead use mean observation parameters to evaluate 
the likelihood on this grid. This precludes the use of FRB 191001, 
which was detected at a lower frequency during commensal obser- 
vations (Bhandari et al. 2020b). The remaining seven FRBs used 
are given in Table 6. 


7 INITIAL RESULTS 
7.1 Calculations 


We use a brute-force approach to find the best-fit values of 
Emax, Y, N, Hhost: Chost and evaluate Eq. 25 over a multi- 
demensional cube of parameter values. The resulting likelihood 
dependence on each is calculated for single (pairs of) parameters by 
marginalising over the remaining five (four) parameters. That is, the 
value taken is the maximum likelihood found when the remaining 
five (four) values are varied over their full range. 

In this work, we use a frequentist approach to setting confidence 
limits. Confidence intervals are determined using Wilks’ theorem, 


2(logLmax—logL) ~ Xy (41) 


where y2,, is a chi-square distribution with ndf degrees of free- 
dom, here equal to the number of parameters which have not been 
marginalised over (either one or two throughout this work; Wilks 
1962). The 51 FRBs used in this work should satisfy the large-N 
limit required for Eq. 41 to be valid. 
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Table 3. Observational properties of follow-up observations, for ASKAP Fly’s Eye (FE) and incoherent sum (ICS) modes (Shannon et al. 2018; Bannister 
et al. 2019); the Parkes multibeam (MB) receiver (Keane et al. 2017); and the Greenbank Telescope’s (GBT’s) 820 MHz primary focus and L-band receivers 
(Kumar et al. 2019). From left to right: the telescope and receiver names, the central frequency ¥ and total bandwidth Av, time- and frequency- resolutions dt 


and dy, and nominal sensitivity to a 1 ms duration burst. 


Telescope Nerp  Tovs [hr] 


Mode [MHz] 


ASKAP/FE 20 26,616 1315 
ASKAP/ICS T N/A 1315 
Parkes/Mb 13 3,946 1382 


Av ot ôv Fi 
[MHz] [ms] [MHz] [Jy ms] 
336 1.2565 1 21.9 
336 1.2565 1 4.4 

337.1 0.064 0.39 0.5 


Table 4. Properties of fast radio bursts detected by the Parkes multibeam 
system, and used in this work. Given is the original FRB designation, mea- 
sured total DM and DMjsy estimated by the NE2001 (Cordes & Lazio 
2002) in pe cm, and ratio of measured to threshold SNR. Entries with a 
* indicates that this value is approximate. 


Parkes: total Tops = 164.4 days 


FRB DM  DMism $ Ref. 
110214 168.9 32 1.44 Petroff et al. (2019) 
110220 944.4 36 5.44 
110627 723 48 1.22 
110703 1103.6 33 izg  eehtoneral, (OUTS) 
120127 553.3 33 1.22 
090625 899.55 32 2.8 
121002 1629.18 72 1.6 
130626 952.4 65 2 Champion et al. (2016) 
130628 469.88 52 2.9 
130729 861 32 1.4 
151230 960.4 48 1.7 , 
160102 2596.1 36 jes Siamani GMa 


Parkes: unnormalised observation time 


FRB DM è DMsm s Ref. 
010305 350 44 1.02* Zhang et al. (2020a) 
010312 1187 51 1.1* Zhang et al. (2019) 
010724 375 45 2.3* Lorimer et al. (2007) 
131104 779 71 3.06 Ravi et al. (2015) 
140514 562.7 36 1.6 Petroff et al. (2015) 
150807 266.5 38 5* Ravi et al. (2016) 
180309 263.52 46 41.1 ; 
180311 1570.9 46 q5, eee GOND 


7.2 Degeneracy 


Initial calculations revealed a degeneracy in the fitting parameters 
between Emax, œ, and n. This is shown in Figure 7, which plots 
the variation of the marginalised likelihood over these parameters. 
In each case, there is a broad maximum extending over the full 
range of a (—5 — —0.5). The reason for this degeneracy is that these 
three parameters are strongly related to the high-DM, high-z cut-off 
in the observed FRB distribution. This is restricted by the lack of 
ASKAP/ICS-localised FRBs at high redshifts, and by the number 
of high-DM FRBs detected by ASKAP/FE and Parkes multibeam 
observations. Having a steep spectral index (very negative value of 
œ) provides a mechanism to reduce the expected number of FRBs 
via the k-correction, and is both consistent with, and requires, a 
higher value of Emax. Similarly, it also allows for a rapidly evolving 
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Table 5. As per Table 4, for ASKAP/FE observations. 


ASKAPIFE: Tops = 1274.6 days 


FRB DM DMism s Ref. 


170107 609.5 37 1.68 Bannister et al. (2017) 


170416 523.2 40 1.38 
170428 991.7 40 1.11 
170707 235.2 38.5 1.00 
170712 312.8 35.9 1.34 
170906 390.3 38.9 1.79 
171003 463.2 40.5 1.45 
171004 304.0 38.5 1.15 
171019 460.8 37 2.46 
171020 114.1 38.4 2.05 
171116 618.5 35.8 1.24 
171213 158.6 36.8 2.64 
171216 203.1 37.2 1* 

180110 715.7 38.8 3.75 
180119 402.7 35.6 1.67 
180128.0 441.4 32.3 1.31 
180128.2 495.9 40.5 1.01 
180130 343.5 38.7 1.08 
180131 657.7 39.5 1.45 
180212 167.5 30.5 1.93 


Shannon et al. (2018) 


ASKAP/FE: unnormalised time 


ERB DM DMjm $8 Ref. 


180417 474.8 26.1 1.84 
180515 355.2 32.6 1.27 


Agarwal et al. (2019) 
Bhandari et al. (2019) 


180324 431 64 1.03 


180525 388.1 30.8 2.93 Macquartetal. (2019) 


Table 6. As per Table 4, for ASKAP/ICS observations, with redshift z 
included. 


ASKAP/ICS: unnormalised time 


FRB DM DMjm 8$ z Ref. 


180924 362.4 40.5 2.34 0.3214 
181112 589.0 40.2 2.14 0.4755 


Bannister et al. (2019) 
Prochaska et al. (2019) 


190102 364.5 57.3 1.38 0.291 
190608 339.5 37.2 1.79 0.1178 
190611.2 322.2 57.6 1.03 0.378 
190711 594.6 56.6 2.64 0.522 


Macquart et al. (2020) 


190714 504.7 38.5 1.19 0.209 Heintz et al. (2020) 
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Figure 7. Degeneracy of the marginalised likelihood with œ. Top: depen- 
dence of Emax, bottom: dependence on n. The contours give the 68%, 90%, 
and 95% confidence limits. The step-like behaviour in Emax is due to coarse 
gridding. 


population with redshift (high n); a very large number of high-z 
events is otherwise excluded when a is near zero. This degeneracy 
was also noted by Lu & Piro (2019), who by default use a = —1.5. 

Without further data, or a prior on any of these three parame- 
ters, this degeneracy cannot be broken. All results from the literature 
on Emax and n derive from similar analyses as in this work, albeit 
with simpler methods, and are therefore not independent. However, 
the work of Macquart et al. (2019), who analyse the instrument- 
corrected spectral structure of bursts detected in ASKAP/FE ob- 
servations, provides an independent constraint on this parameter. 
We therefore use a Gaussian prior on «œ, with mean at a = -1.5, 
and oq = 0.3, reflecting their result of a = ao. In doing 
so, we note that the spectral index model may be interpreted as a 
frequency-dependent rate as discussed in Section 3, and that while 
a spectral break below GHz frequencies is expected (Sokolowski 
et al. 2018), the FRBs used in this analysis were all observed at 
GHz frequencies, and thus the results are only sensitive to spectral 
behaviour at — and due to redshift, above — 1 GHz. 


7.3 Comparison of results 


The limits on single FRB population parameters presented in James 
et al. (2021) are obtained with this approach. We observe that while 
our prior on & shifts the preferred values of the other parameters by 
small amounts, it is not a large influence, and primarily acts to limit 
very strong source evolution models with n > 2. Table 7 compares 
these results with a prior on a to those of other authors, as well as 
a brief summary of which effects are included. 


8 RESULTS 2: REDSHIFT AND DISPERSION MEASURE 


The best fitting parameter set (James et al. 2021) allows a com- 
parison between the expected and observed distributions of FRBs 
in DM, z, and SNR space. However, each combination of allowed 
parameters produces a unique map of the z-DM distribution of 
FRBs for each telescope. Rather than present a very large number 
of plots, we use the following approach to identify a finite number 
of reasonable possibilities. 

For each parameter, we take its value at both the upper and 
lower limits of its one-dimensional 90% C.L., and choose the corre- 
sponding values of the other parameters. This results in 12 further 
parameter sets, which are listed in Table 8. For comparison, we also 
include the best-fit parameter set assuming no source evolution. 

For each of these parameter sets — including the best-fit set 
— we generate the expected FRB distribution in DM-z space, and 
compare this to the observed values. Finally, we also consider the 
best-fit set of parameters when setting n = 0, in order to illustrate 
how well this case fits the data. All the results shown in this section 
are calculated using the spectral index interpretation of «œ. 

In interpreting the results in this section, we caution that the 
data to which we compare expectations have been used to deter- 
mine the FRB population parameters, so the two are not indepen- 
dent. However, given that only eight free parameters (the standard 
six, plus two for the FRB width distribution) are used to fit multi- 
ple observables from three different FRB surveys, good agreement 
should not be the result of over-fitting, but rather indicate a genuine 
correspondence between the models used in this work and reality. 


8.1 Observed and predicted distributions 


The predicted DM and z distributions from all fourteen cases — 
best-fit, no source evolution, and twelve sets reflecting parameter 
uncertainty — are shown in Figures 8 and 9 respectively. In the case 
of the DM distribution, data and predictions are shown individually 
from each survey described in Section 6, and stacked together to 
allow a better comparison. For the z distribution, the only available 
data comes from ASKAP/ICS observations — however, predictions 
from each individual survey are also shown. 

The first question to ask is — are the best fits indeed good 
fits? Our best-fitting parameter estimates (James et al. 2021) do not 
necessarily indicate that the modelled DM, z, and SNR distributions 
are good fits to the data — merely that they are the best fits given the 
form of the model used. A by-eye analysis shows that the predictions 
from the best-fit parameter set are indeed a good match to the data. 
However, there appears to be an over-prediction of FRBs at low 
DM and particularly at low redshift. This is a common feature over 
all the parameter sets within the 90% error margins, although the 
degree of peakedness near z = 0 varies greatly. 

We consider four possible explanations for this below. 
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Author Emax a y n Hhost Chost Beam DM|z 7(DM) f(w) 
Lu & Piro (2019) 43.1411 =1.5¢ -0.640.3  0.3*)-94 N N N N N N 
Luo et al. (2020) 416” 04 -0.794916 04 306e 0.17806 Y Y Y Y 
Arcus et al. (2021) 45.14 -2.111, -0.740-35 11 N N N N Y Y 
Caleb et al. (2016) 41.28 0s N/A 0,14 100% 09 Y Y Y Y 
Macquart et al. (2020) N/A N/A N/A N/A este 0. ze N Y N N 

r +0.49 +0.21 +0.11 +0.25 +64 +0.15 

This work 41.844049  LISSHOZIF  —1.16t011 1,774925  145+%4 0.511013 y Y Y Y 


Table 7. Summary of results and methods used by other authors in estimating FRB population parameters. Shown in columns 2-7 are the best-fit parameter 
values and corresponding uncertainties, converted to the units and notation used in this work. For the sake of comparability, all limits have been converted to 1 o 
standard deviations assuming that uncertainties are Gaussian distributed. Columns 8-11 summarise the different effects that were (Y) or were not (N) included 
in the modelling: accounting for instrumental beamshape, a distribution in dispersion measure for a given redshift, a DM-dependent telescope threshold, and a 
distribution of burst widths f (w) and its corresponding effects on sensitivity. Note that a ‘Y’ does not necessarily mean a comparable treatment to this work. “ 
Parameter held fixed by the authors. ” indicates the use of a Schechter function, with exponential decay after Emax. © The host galaxy contribution is modelled 
according to the “ALGs(YMW 16)” model from Luo et al. (2018), with approximate mean ~ 0.8 and standard deviation 0.2 in logjọ DM, to which is added a 
‘local’ contribution with a uniform distribution from 0-50, for a total mean of 30 and quadrature-added deviation 15, i.e. 0.17 in log)g. d The (1+z) scaling 
used in this work is converted to n via n = 2.78 from Eq. 13. © These values are not explicitly quoted, and were approximately derived from plots shown in 
the text. f This constraint primarily arises from the prior of Macquart et al. (2019). A log-normal luminosity function was used; here we quote the mean of 
10°8-? erg plus three standard deviations of 1 in logy. 


Table 8. Parameter sets used in this work: the best-fit values found in our companion paper (James et al. 2021); sets allowing the parameter in the left-most 
column to take its minimum and maximum value within 90% C.L; the best-fitting parameter set when forcing n = 0, i.e. no source evolution; and a set when 
setting Emin = 1038-5 erg, but otherwise identical to the best fit. 


Set logig Emin 10810 Emax Qa y n Hx Ox 
Best fit 30 41.84 -1.55 -1.16 1.77 2.16 0.51 
logj9 Emax 30.0 41.6 -1.5 -0.8 1.0 2.0 0.9 
30.0 42.51 -1.5 -1.25 1.91 2.0 0.45 

a 30.0 41.9 -1.88 -1.15 1.9 2.25 0.5 

30.0 41.64 -1.2 -1.12 1.4 2.25 0.5 

y 30.0 42.08 -1.5 -1.34 2.08 2.25 0.5 

30.0 41.8 -1.5 -0.96 1.52 2.0 0.6 

n 30.0 41.8 -1.5 -1.1 1.11 2.25 0.5 

30.0 41.88 -1.75 -1.2 2.28 2.15 0.5 

Hhost 30.0 42.18 -1.5 -1.1 1.78 1.77 0.59 
30.0 41.8 -1.5 -1.2 1.67 2.41 0.56 

Chost 30.0 42.16 -1.5 -1.2 1.8 2.08 0.36 
30.0 42.0 -1.5 -1.1 1.6 20 0.81 

n=0 30. 41.6 -1.25 -1.1 0. 2.5 0.6 
Eri 38.5 41.84 -1.55 -1.16 1.77 2.16 0.51 


8.2 Random fluctuations for the DM and z distributions respectively. We therefore conclude 
that the apparent deficit of FRBs at low DM and redshift compared 
with predictions is consistent within 1 ø statistical fluctuations of 
expectations. Nonetheless, we proceed with further analysis since 
the presence or otherwise of a minimum energy, or effects due to a 


large fraction of the population being repeaters, is of great interest. 


The observed deficit of low-DM and low-z FRBs may simply be 
a statistical fluctuation. To evaluate this, we perform Kolmogorov- 
Smirnov tests on the ASKAP ICS z distribution and total DM distri- 
bution (Kolmogorov 1933; Smirnov 1948). Predicted and measured 
cumulative distributions of both DM and z are shown in Figure 10 
— the KS-statistic is the maximum absolute difference between the 


be ; as 8.3 Evidence for a minimum energy Emi 
two curves. To evaluate the significance of this statistic, we take the BY “min 


best-fit curves as the truth, and randomly generate 10,000 samples 
from each. Histogramming the results produces the expected distri- 
butions of the KS statistic under the null hypothesis that the best-fit 
prediction is true. Comparing this distribution to the observed value 
of the KS statistic in Figure 11 shows that our observed values of 
DM and z are completely consistent with predictions, with larger 
values of the KS statistics observed in 17% and 41% of cases for 
the DM and z distributions respectively. Performing a similar anal- 
ysis using the n = 0 distribution gives p-values of 22% and 2% 
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In our standard modelling, we have set the minimum FRB energy 
E min to an extremely low value of 1030 erg — well below the char- 
acteristic energies of observed FRBs, effectively making it zero. 
This is because values of y > —1.5 render FRB observations pri- 
marily sensitive to Emax (Macquart & Ekers 2018b), while bursts 
from FRB 121102 have been observed at much lower energies than 
are likely to be probed by Parkes and ASKAP observations (e.g. 
Law et al. 2017). 

However, a clear possible explanation for the apparent deficit 
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Figure 8. Observed (histogram) and predicted (lines) distributions of DM for all surveys considered (upper left), and the individual Parkes multibeam (upper 
right), ASKAP Fly’s Eye (lower left), and ASKAP ICS (lower right) surveys. Predictions show the best-fit over the entire parameter space, when constrained 
to no source evolution, when allowing Emin to vary, and when varying parameters within their 90% C.L. 


of FRBs at low DM-z is a minimum FRB energy. Without such a 
cut-off, as telescopes probe ever lower values of E in the nearby 
Universe, the predicted number of FRB observations per comoving 
volume will increase without limit when y < —1.If-1 > y > -1.5, 
the reducing volume will somewhat compensate, and the total rate 
will remain finite. This gives rise to the sharp increase in the best-fit 
expected redshift distributions near z = 0 in Figure 9, although such 
a peak may not be present within 90% error margins. 


To investigate the effect of Emin, we fix the best-fit param- 
eters, and vary Emin. The evolution of the likelihoods is shown 
in Figure 12. Interestingly, the total likelihood decreases with in- 
creasing Emin, With Emin = 1038-5 erg the 90% C.L. upper limit. 
Why? As expected from Figures 8 and 9, the p(z, DM) contribution 
increases with Emin, peaking near 1038-5 erg. This peak is a combi- 
nation of the ASKAP/ICS observations strongly favouring a large 
Emin © 1040-5 erg, and the Parkes multibeam observations strongly 
favouring Emin < 103° erg, with the ASKAP/FE observations being 
relatively neutral until Emin > 10% erg. 


The reason why the Parkes observations are strongly against 
high values of Emin is that sensitive telescopes with small fields 
of view are highly unlikely to observe low-DM bursts, unless the 
small volume of the near Universe corresponding to such a low 
DM is populated by many (necessarily low-energy) FRBs. The 
lowest-DM burst detected by Parkes during the included surveys, 
FRB 110214, had a DM of 168.5 pc cm7? (Petroff et al. 2019). Set- 


ting Emin = 1038-5 implies reduced FRB rates for redshifts closer 


than z = 0.194, assuming a limiting fluence of 0.3 Jy ms, with those 
bursts that are detected likely to have a SNR significantly greater 
than threshold. A redshift of z = 0.194 implies a most likely dis- 
persion measure of approximately 281 pc cm3, being composed 
of a cosmological contribution of DMeosmic = 83 pc cm™?, lo- 
cal contribution of 82pccm~?, and our model best-fit value of 
Hnost/(1 +z) = 116pe cm~>. The observation of FRBs by Parkes 
with DMs below this value therefore disfavour a significant mini- 
mum energy Emin- 


This is illustrated in Figure 13, where we show model pre- 
dictions of ps, and the observed values of s, for each Parkes- 
detected FRB. This is done for the best-fit model, and when using 
Emin = 1038-5 erg. While the Emin model predicts that the high 
fluence of FRB 180309 with a DM of ~263 pcem™3 and s = 41 
is slightly more likely, it predicts that FRB 110214, with a DM 
of ~169 pc cm7>, is much less likely to have its observed value of 
s = 1.44. 


In this work, we have not included Emin as a global minimi- 
sation parameter due to computational constraints. Might there be 
some other combination of parameters for which a significant Emin 
is found? To investigate this, we have repeated the Emin optimisa- 
tion for all parameter sets in Table 8. Only for the parameter set 
minimising y do we find a significant value of Emin to be preferred, 
with a best-fit value of 1038-0 erg, and 90% upper limit of 1038-8 erg. 
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Figure 9. Predicted distributions (lines) of z for Parkes multibeam (top) and 
ASKAP Fly’s Eye (middle) surveys, with the ASKAP ICS results (bottom) 
also showing the observed values (histogram). Predictions show the best-fit 
over the entire parameter space, when constrained to no source evolution, 
when setting Emin = 10°° erg, and when varying parameters within their 
90% C.L. 


This makes sense, since otherwise a steep energy function would 
over-predict the number of near-Universe FRBs. The resulting gain 
in likelihood acts to slightly weaken our confidence in the lower 
limits on y, e.g. the 90% lower limit shifts to 81% confidence. 

Our findings clearly do not indicate that there is no minimum 
FRB energy. Within the confines of our power-law model, the data 
used appear insensitive to Emin < 1037 erg, and we can only rule 
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Figure 10. Observed and predicted cumulative distribution functions for the 
DM from all considered surveys (above) and z from the ASKAP ICS survey 
(below). Predictions show the best-fit over the entire parameter space, when 
constrained to no source evolution, and when allowing Emin to vary. 


out Emin > 1038-5 erg at 90% C.L. — but it does run counter 
to the findings of Luo et al. (2020). These authors find a most 
likely minimum luminosity of 104? erg s~!, which is approximately 
1039 erg assuming a standard 1 ms burst, although they conclude 
rather that this finding is due to the limit of their sample. 


8.4 Influence of repetition 


In this study, we have ignored repeating FRBs on the basis that none 
of the FRBs detected by ASKAP and Parkes have been observed 
more than once — even though some are known to repeat (Kumar 
et al. 2019; Patek & Chime/Frb Collaboration 2019; Kumar et al. 
2020), the surveys were not sensitive enough to probe this. How- 
ever, the method of Section 5 covers the entirety of z-DM space, 
regardless of whether or not an FRB has been detected at that point. 
Therefore, for any hypothesised repeating FRB, there will always 
be a sufficiently nearby volume of the Universe where any survey 
would be expected to detect more than one burst if the FRB landed 
in its field of view. If such a repeating FRB happens to be located 
in that volume, the observed burst rate will be greater than expected 
— but if one does not, it will be less. 

This effect is analysed in the context of Canadian Hydrogen 
Intensity Mapping Experiment (CHIME) observations by Gardenier 
et al. (2020), who note that the DM distribution of repeating FRBs 
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Figure 11. Observed value (red line) and simulated distribution (blue his- 
togram) of the KS statistics for the DM (top) and redshift (bottom) distribu- 
tions under the best-fit hypothesis. Shown are the corresponding p-values to 
reject this hypothesis. 
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Figure 12. Evolution of likelihoods as a function of minimum energy Emin, 
relative to the value at Emin = 10% erg. Shown is the total likelihood; the 
likelihood over all surveys, split into contributions from the z-DM distribu- 
tion p(z,DM), the number of events p(N), and the signal-to-noise probability 
p(SNR). The p(z,DM) contribution is further split into components from the 
ASKAP/FE, ASKAP/ICS, and Parkes/Multibeam surveys. 
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Figure 13. Likelihoods of observing FRBs with a given s, ps(s), as a 
function of their dispersion measure, DM, for the Parkes multibeam obser- 
vations, weighted by s for clarity. This is calculated for best-fit parameters 
using Emin = 103° erg (top), and Emin = 1038-5 erg (bottom). 


should be lower than that of repeaters observed only once. While 
this has not yet been observed (CHIME/FRB Collaboration et al. 
2019a,b; Fonseca et al. 2020a,b), this may be due to a very broad 
distribution of intrinsic repetition rates (James et al. 2020b). 

As discussed by James (2019), the ASKAP/FE observations 
have deeply probed some regions of sky with over 50 antenna- 
days spent on individual fields. This has allowed strong limits to 
be placed on FRB repetitions — but also made the observations 
more susceptible to whether or not a strong repeater at z < 0.2 was 
present in these fields. We observe that the low-DM deficit is greatest 
in ASKAP/FE observations, and lowest for Parkes observations, 
consistent with this prediction. 

Until more is known about the population of repeating FRBs, 
we cannot further quantify this effect, except to add that this effect 
is guaranteed to be present to some degree (some FRBs definitely 
do repeat), with its importance increasing if all FRBs are explained 
by strong repeaters, and lessening as a lower fraction of all FRBs 
are due to repeaters, and those repeaters are weak. 
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8.5 Minimum search DM 


Many FRB surveys use a minimum searched DM, either to exclude 
Galactic FRB candidates, or due to local RFI with intrinsic DM of 0 
that nonetheless contaminates searches. For example, Thornton et al. 
(2013) reject bursts with DM< 100 pe cm~?, although in our model 
nearby FRBs could have a DM of as little as DM 80 pc cm3. 
This effect will sometimes exclude FRBs in the very nearby Uni- 
verse, potentially resulting in the observed deficit. Only two bursts 
however — 171020 (Shannon et al. 2018) and 180729.J0558+56 
(CHIME/FRB Collaboration et al. 2019a) — have been detected 
near this limit, and the deficit extends well beyond 100 pc cm3, 
so we consider this an unlikely cause. Analysis of historical Parkes 
data by Zhang et al. (2020b) has found one new FRB at 350 pe cm~? 
(which was too recent to include in this analysis), but none at lower 
DM. 


8.6 Summary: redshift and dispersion measure distributions 


Having established the robustness of these results, we summarise 
the predictions for the redshift distributions from Figure 9. In all 
scenarios, the redshift distribution of the ASKAP/FE detections lie 
in the range zmax < 0.7, and most fits find zmax < 0.5. Over all 
parameter sets, between 32% and 49% of ASKAP/FE bursts should 
originate from within z < 0.1, confirming that these bursts are ideal 
follow-up targets due to their likely proximity. This suggests that 
the limits set on the repeatability of individual FRBs by James et al. 
(2020a) are significantly stronger than published, since those authors 
conservatively assume maximal redshifts. It also lends additional 
weight to the possible association of FRB 171020 with a galaxy at 
40 Mpc by Mahony et al. (2018). 

The predicted z-distribution of Parkes bursts is significantly 
broader than that of ASKAP/FE observations, although the no- 
evolution scenario still suggests a large fraction (35-60% at 90% 
confidence over all parameter sets) in the z < 0.5 Universe, with 
6-20% (over all parameter sets) being within z < 0.1. 

A key test of our prediction of a large number of near-Universe 
FRBs will be future ASKAP/ICS detections. So-far, all bursts de- 
tected by ASKAP’s ICS mode have been in a limited range of both 
DM and z, which seem to be from the central part of all predicted 
distributions. Our best-fit model predicts that 24% of ASKAP/ICS 
localisations should lie within z < 0.1, with a range over all sets 
of 14-32%. To date, none have been observed — even amongst the 
unpublished ones not included here. 


8.7 Source counts distribution 


The slope of the source counts (“logN—logS”) distribution was one 
of the first FRB observables to be analysed. Adapted from its original 
use in the study of radio galaxies characterised by their flux S (Ryle 
1968), applied to FRBs, the source counts distribution is the number 
N of observed FRBs as a function of fluence threshold Fip. In an 
infinite Euclidean Universe, the distribution is expected to have a 
form 


Fin \“ 
N(Fin) = C F7 (42) 
0 
with a = —1.5, and C and Fo normalising constants. Studies us- 


ing a variety of methods, with different treatments — or neglect 
of — observational biases in Fip, and using data from telescopes 
with different detection thresholds, found inconsistent values of a 
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Figure 14. Predicted source counts distribution, p(s), multiplied by s!> 
for clarity, for top: ASKAP/lat50, and bottom: Parkes/Mb, according to 
the best-fit parameter set (black solid), when setting Emin = 10% erg (red 
dashed), and when varying parameters within their 90% C.L. (thin grey). 
Also shown are lines of constant power-law slopes to guide the eye. The 
fluctuations above s = 10 are due to finite gridding in z and DM. 


in the range of —0.8 > a > —2.6 (Vedantham et al. 2016; Opper- 
mann et al. 2018; Lawrence et al. 2017; Macquart & Ekers 2018a; 
Bhandari et al. 2018). James et al. (2019b), reviving methods ap- 
plied to radio galaxy studies by Crawford et al. (1970), argue that 
s — defined in Eq. 34 — is a bias-independent measure of the 
source-counts slope, and find a = —1.18 + 0.24 (68% C.L.) for 
Parkes FRBs, and a = —2.2 + 0.47 for ASKAP/FE FRBs, equating 
to a 2.60 tension. This was qualitatively consistent with Macquart 
& Ekers (2018b), who argue that at high values of Fih, the slope 
should be Euclidean (a = —1.5), while the parameters of the FRB 
population will lead to a flattening at lower thresholds. 

Our model can be used to derive the expected distribution of 
s using Eq. 36, by integrating over all values of DM and z, then 
converting this differential distribution to a cumulative distribution. 
The results for ASKAP/FE and Parkes/Mb observations are given 
in Figure 14. 

We immediately see that in all scenarios, ASKAP/FE obser- 
vations, with a higher base threshold of 26 Jy ms, are expected to 
follow a Euclidean (a ~ —1.5) distribution, while near-threshold 
(1 < s < 10) Parkes/Mb observations exhibit a flatter source-counts 
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index near a = —1.3. Different scenarios predict source-counts in- 
dices in the range -1.4 < a < —1.6 above s = 10 for Parkes/Mb 
observations. The existence or otherwise of a minimum energy does 
not significantly affect the distribution. 

Our results in all scenarios are consistent with the findings 
of James et al. (2019b), with the greatest tension — about 1.50 
— being between the source-counts slope that those authors find 
for ASKAP of —2.2 + 0.47 (68% C.L.), and the approximate range 
of —1.4 to —1.6 found here. In particular, we also find that the 
source-counts index for Parkes/Mb is expected to be flatter than for 
ASKAP/lat50, and furthermore, the apparent ‘deficit’ of FRBs with 
low values of s is potentially attributable to the true behaviour of 
the FRB population, rather than statistical fluctuations or a mea- 
surement bias. This also suggests that the lower values of a found 
by previous authors — Vedantham et al. (2016); Oppermann et al. 
(2018); Lawrence et al. (2017) — incorporating data from tele- 
scopes more sensitive than Parkes may have been correct. 


9 THE Z-DM DISTRIBUTION 


We argue in this work that the best representation of the FRB pop- 
ulation observable by a telescope is a two-dimensional function 
of extragalactic (cosmological plus host) dispersion measure, and 
red shift. Necessarily, the observable fraction of this distribution is 
a function of survey parameters, and also local DM contribution, 
which reduces sensitivity as it increases at lower Galactic latitudes. 
Our best-fit z-DM distributions for the three surveys considered are 
plotted in Figure 15. 

We describe the general features of these plots. Most FRBs are 
expected to lie near to, or below, the Macquart relation, being the 
1-1 correspondence of FRB redshfit with DM. This relation, with a 
slope of approximately 100 pc cm3 per 0.1 in redshift, continues up 
to some maximum detectable distance, being z=0.5, 0.9, and 2.4 for 
ASKAP/FE, ASKAP/ICS, and Parkes/Mb respectively. However, 
the majority of FRBs will not be found near the maximum redshift, 
simply because there are very many more FRBs with lower energy, 
and much more of the sky covered at lower sensitivity. This is most 
evident with Parkes/Mb observations, where the most likely half of 
all FRB observations will arise on the Macquart relation with z < 1. 

Off the Macquart relation, there is a significant fraction of 
FRBs expected to be found with higher than expected DMs, due to 
a combination of their host and cosmological contributions. Only 
our ASKAP/ICS (localised) FRB sample has provable examples, be- 
ing FRB 190714 (504.7 pecm™? from z = 0.209) and FRB 190608 
(339.5 pe cm~? from z = 0.1 178), although a third (FRB 191001, 
with 506.92 pc cm~? from z = 0.23) was excluded due to being de- 
tected at a lower frequency. This effect is less pronounced for more- 
sensitive surveys, since excess host contributions are dominated by 
cosmological ones. We emphasise that much of the structure above 
the Macquart relation — i.e. in the high-DM region — is poorly 
constrained, since our adopted log-normal distributions may not 
reflect reality. 

At very high DMs, only near-Universe FRBs are observable, 
since a burst must be observed with very high fluence to overcome 
the detection bias against high DM. The upper bound of the 99% re- 
gion (dashed lines) slopes backward, against the Macquart relation, 
since more low-energy FRBs with large excess DM are predicted 
than high-energy FRBs lying on the Macquart relation. We do not 
expect our quantitative estimates in this region to be accurate until it 
is directly probed with localised FRBs. Since the cause of this effect 
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Figure 15. Normalised probability distribution p(z, DMgg) of the redshift 
z and extragalactic dispersion measure DMgg of FRBs detected in ASKAP 
Fly’s Eye (top), ICS (middle) and Parkes multibeam (bottom) surveys. Con- 
tours indicate the intervals in which 50% (dotted), 90% (dot-dash), and 99% 
(dashed) of all FRBs would be detected. Red dots indicate the values of 
DMgg and host galaxy redshift z for ASKAP ICS localized FRBs. 
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is a well-understood observational bias however, it will clearly be 
present to some degree. 


9.1 The Macquart Relation 


The ‘Macquart Relation’ is the general one-to-one-ness of the rela- 
tionship between redshift and DM of FRBs. It is predicted from the 
distribution of baryonic matter in the Universe (Inoue 2004), and 
first evidence for it was given in Shannon et al. (2018), by com- 
paring the DM distributions of ASKAP/FE and Parkes/Mb popu- 
lations, where the higher sensitivity of Parkes allowed it to probe 
more-distant FRBs with higher DMs. The relation was first observed 
directly by Macquart et al. (2020), who showed that the redshifts 
of localized FRBs were consistent with the baryonic content of the 
Universe. 

For the purpose of comparing survey results, we propose a 
useful distinction: the ‘weak’ Macquart relation (which might more 
accurately be titled the ‘Shannon’ relation), where telescopes with 
higher sensitivity on average observe more-distant FRBs with higher 
DM; and the ‘strong’ version (or true Macquart relation), where the 
DMs of FRBs within a survey are a good proxy for their redshift. 
Several authors use a 1-1 z-DM relation in performing estimates of 
the FRB population from the DMs of un-localised FRBs (Shannon 
et al. 2018; Deng et al. 2019; Lu & Piro 2019; James et al. 2020a). 

Clearly, the weak version of the relation is well-established, 
and an obvious consequence of the cosmological nature of FRBs. 
Here, we test the strong version, which can be readily tested by 
calculating 


p(z,DM) 
p(DM) 


for a given survey. This is shown for each survey in Figure 16. 

In each survey, the Macquart relation applies up to a maxi- 
mum redshift zmax, beyond which it reverses. The reversal can be 
intuitively understood by realising that at the maximum redshift at 
which an FRB can be detected, FRBs with any excess of DM can not 
be detected, due to DM smearing reducing sensitivity. Therefore, 
the only way to detect FRBs with a DM lying above that expected 
from a burst originating at Zmax is to have the burst originate at 
a nearer redshift. As noted in Section 3 and Fig. 2, the increased 
number of FRBs in the local Universe is related to the cumulative 
luminosity index y = —1.16, with more negative values leading to 
more nearby high-DM events. 

The reversal of the Macquart relation has several practical 
consequences. Firstly: for surveys with a large sample of FRBs, 
the burst with the greatest DM will not be the most distant. An 
excellent example of this phenomenon is FRB 170428 (ASKAP; 
Shannon et al. 2018), which is most likely to originate below z = 0.3, 
rather than the value of z ~ 1 expected from the Macquart relation. 
FRB 160102, observed by Parkes, is another likely candidate. The 
implication is that works using a 1-1 DM-~ relation will vastly 
over-estimate the maximum FRB energy, since they will necessarily 
attribute a large distance and therefore high energy to the highest 
DM event, which may in fact be quite local. 

A key test of this reversal would be the detection of an FRB with 
DMz 1000 pe cm? by ASKAP in ICS mode, and its subsequent 
localisation to a redshift z < 0.5. Again, we note that this reversal of 
the Macquart relation will always be present to some extent, since it 
is fundamentally due to observational effects which are known and 
understood — it is merely the extent of this effect which is currently 
uncertain. 

Finally, we note that the existence of FRBs with very high DMs 


P(z|IDM) = (43) 


MNRAS 000, 1-23 (2020) 


FRB z—DM distribution 19 


0.75 1.0 
20 
lo 
0.6 Median oe 
= 
Macquart a 
N 
0.45 0.6 ~ 
£ 
i g 
= 
0.3 0.4 = 
N 
Q 
0.15 0.2 
0.0 a 0.0 
0. 200. 400. 600. 800. 1000. 
DMeG 
1.0 1.0 
20 
a lo 
0.8 F A Median 0.8 = 
aba a 
= N 
0.6 T =f | 0.6 E 
= JA é 
ag z 
0.4 0.4 A 
N 
Q 
0.2 0.2 
0.0 2 — r 0.0 
0 400 800. 1200. 1600. 2000. 
DMeG 
1.0 
lo 
Median KEES na 
N = 
eect a 
E N 
0.6 E 
h Š 
= 
0.4 a 
N 
Q 
0.2 
0.0 


“0. 800. 


1600. 2400. 3200. 4000. 
DMecG 


Figure 16. Probability p(z|DM) of an FRB originating from redshift z 
given it has been observed with the specified DM in the ASKAP/FE (top), 
ASKAP/ICS (middle), and Parkes/Mb (bottom) surveys. 


has raised the possibility of using FRBs to probe for the signature of 
Helium reionisation (Deng & Zhang 2014; Caleb et al. 2019; Linder 
2020). While this is by no means ruled out, it emphasises that doing 
so will require FRBs to be localised, since simple measures of FRB 
properties as a function of DM will yield a very large scatter in 
redshifts, and hence reduced statistical power. 
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10 CONCLUSION 


We have developed a precise model of FRB observations, including 
observational biases due to the full telescope beamshape, degrada- 
tion in efficiency due to DM, and intrinsic burst width. None of 
these effects are fundamentally new; many others should take credit 
for highlighting their importance, and Luo et al. (2020) should be 
attributed with a first analysis using these techniques. Here, we 
have improved upon this method by using an unbiased data sample, 
adding new observations of localized FRBs, studying the effects of 
source evolution, including the likelihood of the observed signal 
strength, and improving the beam model. We show that ignoring, 
or incorrectly modelling, these factors leads to significant biases in 
the expected redshifts of observable FRBs. 

We have used our approach to model FRB observations with 
ASKAP in both fly’s eye and incoherent sum mode, and Parkes 
multibeam observations. We have carefully selected our data to 
ensure it is not biased due to under-reporting of observation time, or 
due to large local DM contributions reducing sensitivity. Crucially, 
we have included a sample of localized FRBs from ASKAP for 
which the redshift of the host galaxies is measured. 

The z, DM, and SNR distributions of FRBs predicted by our 
best-fit population estimates, presented in James et al. (2021), are 
tested against observations. We find no evidence for, and some 
evidence against, a lower bound to the FRB energy distribution, 
although we only exclude Emin > 1038-5 erg (90% C.L.). No such 
minimum energy is expected for the magnetar-origin hypothesis, 
which links the observed extragalactic FRB population to radio 
bursts from magnetar flares in our own Galaxy. 

Our model also allows us to make inferences on the redshifts of 
the un-localized samples of FRBs detected by ASKAP and Parkes. 
We find these to be somewhat lower than expected from the Mac- 
quart relation, and indeed that the highest-DM events are likely not 
the most distant, due to observational effects causing an inversion 
of the Macquart relation, and the relatively steep best-fit value of 
the FRB energy distribution. The ability of this model to place pri- 
ors on the expected redshift of FRBs given their measured DMs 
will also aid FRB localisation efforts, especially for those bursts — 
such as FRB 171020 — which have uncertain, but promising, host 
associations. 

For the first time, we have incorporated the measured signal- 
to-noise ratio s into FRB population modelling, allowing the use of 
this observable to constrain population parameters, and to predict 
the source-counts (“logN-logS”) distribution. In all scenarios, we 
find a steepening of this distribution from Parkes to ASKAP, con- 
sistent with the predictions of Macquart & Ekers (2018b) and the 
observations of James et al. (2019b). 

Anomalies between observations and our model include a lack 
of ASKAP incoherent sum (localised) FRBs within z < 0.1, an 
over (under) prediction of FRBs from Parkes multibeam (ASKAP 
fly’s eye) observations, and an over-prediction of low-DM bursts. 
We note several potential explanations for these discrepencies, al- 
though none of them have a high statistical significance. Future 
observations of localised FRBs with very large excess DMs, and/or 
from the local Universe, would verify our predictions. FRB surveys 
with very sensitive radio telescopes such as the Five-hundred-meter 
Aperture Spherical Telescope (FAST), or repeating previous Parkes 
and ASKAP surveys at different frequencies, would help to further 
constrain FRB population models. In particular, the application of 
this model to the large sample of bursts observed by CHIME would 
be particularly useful, although it would then need to be adapted 
to include repeating FRBs. We also aim to investigate improved 


numerical/computational methods to speed up calculations to allow 
the inclusion of further data. 
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APPENDIX A: NEGLECTED EFFECTS 


Here we discuss effects that affect search sensitivity which have 
been neglected in this work. 


A1 Fine sensitivity effects 


The precise sensitivity of an incoherent FRB search to an FRB 
depends not only on the ‘coarse’ effects of effective pulse width, but 
also the exact alignment of a burst with the time-frequency binning 
of the data compared to the FRB arrival time, and the DMs which 
are searched compared to the exact DM of the FRB. This effect is 
analysed by Keane & Petroff (201 5a) for early FRB search methods, 
finding fluctuations of up to +15%, although in current FRB search 
methods fluctuations are essentially zero. This was also found to be 
the case for an internal investigation into the CRAFT FRB searches 
with ASKAP for Shannon et al. (2018), with fluctuations of +3%. 
Such effects are thus ignored in this work, and likely should be in 
all future works. 

The presence of radio-frequency interference (RFI) during ob- 
servations can result in a loss of effective bandwidth — and hence 
sensitivity — via vetoed frequency channels; loss of effective ob- 
serving time if the RFI results in FRBs being completely unob- 
servable; and reduced sensitivity to FRBs in a certain parameter 
space (particularly low-DM or high-width FRBs). The sensitivity 
of ASKAP to such effects has been studied using pulsar calibration 
observations by James et al. (2019a), finding a typical 10% loss of 
effective observation time and ~10% fluctuation in sensitivity. 

Further details of search pipelines can effect burst sensitiv- 
ity. The zero-DM subtraction method — i.e. subtracting the mean 
detected power prior to dedispersion — will reduce sensitivity to 
low-DM events and very bright bursts, although typically only if the 
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total dispersion sweep across the detection bandwidth is low. For 
instance, we have found a few percent bias in the estimated SNR 
for ASKAP searches when the expected SNR is greater than one 
hundred. An analysis of the HEIMDALL software — used in most 
Parkes FRB searches — using real-time injected FRBs by Gupta 
et al. (2021) finds a 20% reduction in SNR compared to expecta- 
tions. This mostly affected bursts with widths above 20 ms however, 
and we find in Section 4 our estimates of the DM-z distribution to 
be insensitive to the population of very wide bursts. For ASKAP, the 
expected SNR of pulsar bursts was found by Shannon et al. (2018) 
to be within a few percent of that found by the Fast Real-time En- 
gine for Dedispersing Amplitudes (FREDDA) algorithm used by 
the CRAFT collaboration. 

We also note that in this work we use SNR as the true SNR of 
the FRB, when by definition this has a lo error of unity. Murdoch 
et al. (1973) find that for SNR above six, the measured SNR can 
be approximated as the true SNR for statistical purposes, and so we 
have here ignored this effect, rather than marginalising over it. 

Our beamshape model for Parkes uses the mean observation 
frequency, while for ASKAP it is uniformly weighted over all fre- 
quencies. A preference for e.g. low-frequency bursts would result 
in a slight increase in the survey effective area for both instruments, 
with only the difference in increase being relevant when estimating 
population parameters other than the absolute FRB rate. We con- 
sider this a third-order effect for the frequencies and bandwidths 
considered here. 

Finally, we have not discussed errors in Fo. The absolute scale 
of the mean ASKAP threshold was calculated by James et al. 
(2019a) by referencing observations to Parkes, and using Hydra 
A as an absolute flux calibrator. Thus we expect errors in flux cali- 
bration for Parkes and ASKAP to be linked, and cancel to first order. 
However, individual antenna sensitivity for ASKAP was found to 
vary by ~ +5%, and this is similar to the level of uncertainty 
found when performing calibration observations for Bannister et al. 
(2019). Perhaps the largest source of uncertainty in Fg comes from 
the quoted threshold of 0.5 Jy ms for Parkes, which has an inher- 
ent rounding uncertainty of +10%, and thus a corresponding rate 
uncertainty of +15% for a cumulative source counts index of —1.5 
(see Section 8.7). 


A2 Validity of the measurement of a 


We take as given the measurement of the spectral index by Mac- 
quart et al. (2019). However, as with all other observables used in 
this work, it is valid to ask: why do we not fit the observed dis- 
tribution, taking selection effects into account? This is in-general 
important, since telescope beam widths are broader at low frequen- 
cies, and hence more low-frequency FRBs will be observed. While 
these authors corrected observed FRBs for frequency-dependent 
sensitivity effects, they did not use such effects to model the FRBs 
they did not observe, and the result is thus susceptible to biases due 
to selection effects. 

We first note that in the spectral index interpretation of a, there 
is no such effect — all FRBs are broadband, and every measured 
FRB tends to have the same spectral index. In the rate interpretation 
however, where every FRB is narrow-band, each burst is observed 
according to the sensitivity at its particular frequency. For exper- 
iments with widely spaced beams, the observed sky area scales 
as f -2 so that far more low-frequency bursts would be observed 
even for a frequency-independent rate. In such a case, the observed 
frequency-dependence rate should be used as a fitted parameter. 

The measurements of Macquart et al. (2019) however were 


made with ASKAP, with closely spaced beams overlapping at the 
half-power points. Making a toy calculation, using a Euclidean 
source-counts such that rate is proportional to beam power!-5, we 
find that in the narrow-bandwidth limit, one would find an FRB rate 
proportional to f -1 and hence a = 1. We emphasise that ASKAP 
FRBs are absolutely not in this narrow-bandwidth limit — while 
they do show narrow-band features, these tend to be spread across 
the spectrum (Shannon et al. 2018). We thus expect the bias due to 
this effect on the measurement of œ by Macquart et al. (2019) to be 
very small, and since «œ is a small influence on estimates of other 
parameters in the study, do not consider this further. It would be 
important for CHIME, for instance. 


A3 Pointing 


Here we treat a given FRB survey as having a constant local DM 
component, DMjoca =DMmw+DMhpost. This is because a greater 
local DM reduces sensitivity to extragalatic FRBs due to DM smear- 
ing. 

When this effect becomes significant, this requires extending 
the integral in Eq. 29 to 


<Ni> = Dita [aoa E 
k 


dDMgcp(DMgglz) 


dQdz 
J wow) f awp pE > Ew) (Al) 
Ep ~ Ep(B,w,z, DM = DMgg + DMiocal,k) (A2) 


where 7;, is the total time spent observing at DMjocai = DMjocal,k- 


A4 Scattering 


In Section 4, we combine the intrinsic burst width wint and scatter- 
broadening width wscat into the incident width wine. However, while 
Wint Can reasonably be assumed to be independent of other FRB 
properties, it is plausible that wine will be correlated with both 
DM and/or z, and it will certainly be frequency-dependent. This is 
due to interstellar scattering, which is implicitly included through 
modelling of observed FRB widths, but is not explicitly accounted 
for as per e.g. Caleb et al. (2016). We largely avoid this question in 
this work by choosing surveys of similar frequency — however, we 
consider evidence for such a correlation here. 

Most FRB searches have time resolutions in the range of 
100 us-1 ms, and resolving scattering tails from the intrinsic burst 
structure is difficult: Qiu et al. (2020) are able to do this for only six 
of the 33 FRBs observed at the typically 1.26 ms real-time resolution 
of ASKAP FRB searches. The broader time-frequency structure ex- 
hibited by many repeating FRBs (Fonseca et al. 2020a,b) however 
can often be resolved at these resolutions (e.g. Hessels et al. 2019), 
while observations down to 10 ys have revealed a 1 GHz scattering 
time of only 24 us in FRB 121102, which is obtained indirectly 
from the measured scintillation bandwidth (Michilli et al. 2018). 

The most reliable way to resolve these two contributions is to 
use time-frequency data at the Nyquist resolution. ASKAP (Cho 
et al. 2020; Day et al. 2020) and UTMOST (Farah et al. 2018, 
2019) have analysed such data for six FRBs each. Bursts were found 
to have strong sub-burst structures down to 50 us, with intrinsic 
widths as low as 10 ms. Scattering was conclusively measured for 
a total of four ASKAP FRBs, being in the range 40 ws—3.3 ms at 
1.27 GHz, while for 181112, it was at most 20 us. UTMOST FRBs 
had 835 MHz scattering times of 4 us-30 ms, with one upper limit 
at 0.2 ms. 
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Clearly, there is a broad distribution of FRB intrinsic widths 
and scattering times. There are general effects that this distribution 
can have on the DM-z distribution of FRBs. 

The most complicated potential interaction of scattering and 
the intrinsic width is one that is dependent on the exact position of 
the FRB in DM-z space. 

If indeed FRBs do arise from two source populations with dif- 
ferent cosmological evolutions or different host galaxy properties, 
or elsewise one population of object with properties that age on 
cosmological timescales, the intrinsic width distribution may have 
some DM-z dependence. This possibility should not be ignored. 
However, the most likely redshift dependence arises from the scat- 
tering term. 

Theoretical studies have examined expectations for scattering 
of FRBs during cosmic propagation, with effects attributed to the 
intergalactic medium (IGM), intra-cluster medium (ICM), and the 
halos and interstellar medium (ISM) of intersected galaxies. 

The general form of the redshift dependence is analysed by 
Macquart & Koay (2013), finding that scattering due to the IGM 
will increase as (1 + z)? to z < 1, and as (1+ z)0-2-0.5 for z = 1. 
This is in contrast to the contribution from hosts, which scales as 
(1 + z)~3. Results on the absolute magnitude of the scattering de- 
pend on the assumed minimum and maximum length scales of the 
turbulence. Macquart & Koay (2013) argue that for realistic turbu- 
lence parameters, the total amount of scattering from the IGM and 
ICM is expected to be low, at < 1 and < 5 ms at 300 MHz, respec- 
tively. Zhu et al. (2018) simulate FRBs propagating in a clumpy 
IGM. Examining the dependence of the mean scattering time 7 on 
DM, they find 7 ~ DM!~? for voids, clusters, and filaments over a 
range of simulation parameters. They also conclude that unreason- 
ably high turbulence scales would be required to achieve scattering 
values comparable to that observed in FRBs. 

In some cases, FRBs localisations have been sufficient to iden- 
tify intersections of the line-of-sight with Galactic halos, placing 
upper limits on the degree of scattering caused by such intersec- 
tions (Prochaska et al. 2019; Connor et al. 2020). 

Qiu et al. (2020) find no evidence however for a DM- 
dependence of scattering. Both measurements and expectation 
therefore suggest that the observed distribution of scattering arises 
from the host galaxies. We do not however attempt to model this in 
this work, since the (1 + z)? reduction in width at high redshift will 
be in any case insignificant against the DM smearing effect. A full 
model of scattering will however become important when including 
observations over a wide frequency range. 


A5 Influence of specific parameters 


This work uses measured FRB detection numbers, dispersion mea- 
sures, redshifts, and strengths from three different FRB searches. It 
is useful to probe the influence of each of these on our overall result. 

To do so, we show in Figure Al our confidence limits for 
each parameter using our prior on a when removing the likelihood 
corresponding to py 27, ppmz (eq.31), ps (eq.36), and when 
removing ASKAP/FE, ASKAP/ICE, and Parkes/Mb observations. 
We do this under the spectral index interpretation of œ — results 
are similar under the rate interpretation. 

Unsurprisingly, the most dramatic effect occurs when removing 
information on measured DM and z. In this limit, we do not find 
lower bounds on Emax in the range log jg Emax [erg] > 40, and limits 
on other parameters are very broad. This is not a very interesting, 
informative, or unexpected result, and we ignore the “No ppmz” in 
the following discussion. However, the fact that some constraints are 
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nonetheless derived is evidence for the source-counts distribution 
of FRBs containing useful information (Macquart & Ekers 2018b). 

For all other parameters, we obtain a variety of modified lim- 
its. Clearly, the ASKAP/ICS sample is most constraining to Emax, 
while without the ASKAP/FE sample, the constraints are narrower. 
This is natural, since the ASKAP/FE sample probes the most intrin- 
sically luminous bursts, without which a fit for a narrow population 
distribution is possible. 

The slope of the cumulative luminosity function y is con- 
strained to by flat by pn, without which a steeper (lower) value 
would be obtained, while both the Parkes/Mb sample, and Psnr, pre- 
vent flat (higher) values. This is a very interesting result. As noted 
in our companion paper, our best fit under-predicts the number of 
ASKAP/FE FRBs (12.9 vs. 20 in 1274.6 days), and over-predicts 
those for Parkes/Mb (19.2 vs 12 in 164.4 days). One would expect 
that removing either Parkes/Mb or pn would have a similar effect by 
removing this tension. However, removing Parkes/Mb also removes 
the Parkes DM distribution. This then suggests that the Parkes DM 
distribution is indicative of a steep source-counts spectrum, whereas 
the number of FRBs detected by Parkes is indicative of a flat spec- 
trum. We do not have a full solution to this quandry, although it 
is clearly related to the discussion in Section 8. We do not elab- 
orate further, since this discrepancy is of only marginal statistical 
significance. 

Perhaps our most significant result — strong evidence for 
source evolution with redshift, i.e. n > 0 — is most strongly dis- 
favoured by pn, such that when this is removed, we prefer a strongly 
evolving population. This, predictably, is related to the low number 
of bursts observed by Parkes/Mb relative to ASKAP/FE compared 
to expectations discussed above, which alone argues against source 
evolution. Removing any of the three FRB surveys, or Psnr, sim- 
ply results in mildly less evidence for n > 0 — there is no single 
result strongly favouring n > 0. We thus conclude that our conclu- 
sion on n robust, especially considering uncertainties discussed in 
Appendix A1. 

Our host galaxy parameters, Hhost and Chost — recall these are 
in log DM space — are most affected by the highest- and lowest- 
sensitivity surveys, ASKAP/FE and Parkes/Mb. Without these ob- 
servations, very low values of Hhost (and hence very large values of 
Chost) become possible. This validates the use of DM distributions 
only when studying these parameters. Neither psnr nor pn have a 
significant effect. 


This paper has been typeset from a TEX/IATRX file prepared by the author. 
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Figure A1. Resulting 68%, 90%, and 95% confidence limits (respectively: blue, upper; green, middle; red, lower) on individual parameters under the spectral 
index interpretation of œ, using in descending order: all available information, removing information on the number of detected FRBs, their DMs and redshifts, 
their measured signal-to-noise ratios, and data from ASKAP/FE, ASKAP/ICS, and Parkes/MB, respectively. 
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